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General Features of Supersymmetric Signals at the ILC: Solving the LHC Inverse 

Problem 
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We present the first detailed, large-scale study of the Minimal Supersymmetric Standard Model 
(MSSM) at a y's = 500 GeV International Linear Collider, including full Standard Model back- 
grounds and detector simulation. We investigate 242 points in the MSSM parameter space, which 
we term models, that have been shown by Arkani-Hamed et al. to be difficult to study at the LHC. 
In fact, these points in MSSM parameter space correspond to 162 pairs of models which give in- 
distinguishable signatures at the LHC, giving rise to the so-called LHC Inverse Problem. We first 
determine whether the production of the various SUSY particles is visible above the Standard Model 
background for each of these parameter space points, and then make a detailed comparison of their 
various signatures. Assuming an integrated luminosity of 500 fb~^, we find that only 82 out of 242 
models lead to visible signatures of some kind with a significance > 5 and that only 57(63) out of 
the 162 model pairs are distinguishable at 5(3)a. Our analysis includes PYTHIA and CompHEP 
SUSY signal generation, full matrix element SM backgrounds for all 2 2 , 2 — > 4, and 2 — > 6 
processes, ISR and beamstrahlung generated via WHIZARD/GuineaPig, and employs the fast SiD 
detector simulation org.lcsim. 
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INTRODUCTION 

The LHC is scheduled to begin operations within a 
year and is expected to change the landscape of parti- 
cle physics. While the Standard Model (SM) does an 
excellent job describing all strong interaction and elec- 
troweak data to date |^, there are many reasons to 
be dissatisfied with the SM. Chief among them are is- 
sues related to electroweak symmetry breaking. As is 
by now well-known, the SM with a single Higgs doublet 
that is responsible for generating the masses of both the 
electroweak gauge bosons and fermions encounters diffi- 
culties associated with stability, fine-tuning, and natural- 
ness. Addressing these issues necessitates the existence 
of new physics at the Terascale. To this end, numerous 
creative candidate theories that go beyond the SM have 
been proposed and many yield characteristic signatures 
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at the LHC. When the ATLAS and CMS detectors start 
taking data at the LHC, they will explore this new terri- 
tory. They will then begin the process of identifying the 
nature of physics at the Terascale and of determining how 
it fits into a broader theoretical structure. 

Of the several proposed extensions of the SM that re- 
solve the issues mentioned above, the most celebrated 
is Supersymmetry (SUSY) ||]. Our working hypothesis 
in this paper is that SUSY has been discovered at the 
LHC, i.e., that new particles have been observed and it 
has been determined that they arise from Supersymme- 
try. Identifying new physics as Supersymmetry is in it- 
self a daunting task, and we would be lucky to be in such 
a situation! However, even in this optimistic scenario, 
much work would be left to be done as SUSY is a very 
broad framework. We would want to know which version 
of SUSY nature has realized and for this we would need 
to map the LHC observables to the fundamental parame- 
ters in the weak scale SUSY Lagrangian. A question that 
would arise is whether this Lagrangian can be uniquely 
reconstructed given the full set of LHC measurements. 
This issue has been recently quantified in some detail by 
the important work of Arkani-Hamed, Kane, Thaler and 
Wang (AKTW) Q|, which demonstrates what has come 
to be known as the LHC Inverse Problem. AKTW found 



that even in the simplest reahzation of Supersymmetry, 
the Minimal Supersymmetric Standard Model (MSSM), 
such a unique mapping does not take place given the LHC 
observables alone and that many points in the MSSM pa- 
rameter space cannot be distinguished from each other. 
Here, we extend their study and examine whether data 
from the proposed International Linear Collider (ILC) 
can uniquely perform this inverse mapping and resolve 
the model degeneracies found by AKTW. 

In brief, AKTW considered a restricted Lagrangian pa- 
rameter subspace of the MSSM. They forced all SUSY 
partner masses to lie below 1 TeV (in order to obtain 
a large statistical sample at the LHC), fixed the third 
generation ^-terms to 800 GeV and set the pseudoscalar 
Higgs mass to be 850 GeV. Points in the MSSM param- 
eter space, hereafter referred to as models for brevity, 
were generated at random with the conditions that tan (3 
lies in the range 2-50, squark and gluino Lagrangian 
mass terms lie above 600 GeV, and Lagrangian mass 
terms for the non-strongly interacting particles be greater 
than 100 GeV. 43,026 models were generated in this 15- 
dimensional parameter space under the assumption that 
all parameter ranges were uniformly distributed, i.e., flat 
priors were employed. No further constraints, such as the 
LEP lower bound on the Higgs mass ^ or consistency 
with the relic density of the universe were applied. For 
each model, PYTHIA [| was used to calculate the result- 
ing physical SUSY spectrum and to generate 10 fb^^ of 
SUSY 'data' at the LHC, including all decays and hadron 
showering effects. This 'data' was then piped through the 
PGS fast detector simulation |^ to mimic the effects of 
the ATLAS or CMS detectors. From this 'data,' AKTW 
constructed a very large number of observables associ- 
ated with the production and decay of the SUSY part- 
ners. No SM backgrounds were included in their 'data' 
sample. AKTW then observed that a given set of values 
for these observables along with their associated errors, 
i.e., a fixed region in LHC signature space, corresponded 
to several distinct regions in the 15-dimensional MSSM 
parameter space. This implies that the mapping from 
experimental data to the underlying theory is far from 
unique at the LHG. Normally phenomenologists deter- 
mine the values of observables for a given point in the 
parameter space of some model. Thus, the task of exper- 
imentalists, to determine parameters of the underlying 
theory from the data, can be viewed as an inverse map- 
ping. The analysis of AKTW found that this mapping is 
many-to-one, that is, the underlying parameters cannot 
be uniquely determined. Thus there is an "LHC Inverse 
Problem." Clearly, if one incorporates the existing SM 
backgrounds as well as systematic effects into this kind 
of study, the number of possible models that share in- 
distinguishable signatures will only increase, potentially 
significantly. The LHC Inverse Problem is thus a very 



serious one. It is also important to be reminded that 
AKTW was a general MSSM study in the sense that no 
assumptions were made about the SUSY breaking mecha- 
nism; almost all LHC studies do make such assumptions, 
and therefore study potential signatures only within a 
few restricted scenarios, which may or may not be real- 
ized in Nature. A general MSSM parameter scan gives 
us a much better feeling for the range of possibilities in 
SUSY scenarios. 

However, the fact that an LHC Inverse Problem exists 
is not overly surprising and the real issue we face is how 
to resolve it. In this paper we will begin to address the 
question of whether the models that AKTW found to be 
indistinguishable at the LHC can be resolved by a high 
luminosity e^e~ collider operating at 500 GeV in the cen- 
ter of mass with a polarized initial electron beam, i.e., 
the ILC. Traditional ILC lore indicates this is the case, 
as studies have shown [0, ||], e.g., the mass and couplings 
of any kinematically accessible weakly interacting state 
should be measured at the 1% level or better at the ILC. 
Such precise determinations imply that decay signatures 
and distributions produced by new particles such as the 
SUSY partners will be observed with relative cleanliness 
and be well measured. The LHC Inverse Problem pro- 
vides us with a unique opportunity to test this lore over a 
very wide range of the MSSM parameter space by com- 
paring the signatures of hundreds of models. We will 
show that, as believed, the ILC can generally distinguish 
models, at least in the case of this restricted scenario 
of the MSSM, and we will explore the reasons why it 
fails when it does. We will find that some SUSY mea- 
surements are more difficult to obtain than previously 
thought, and we will identify some problematic areas of 
the MSSM parameter space which require further study. 
We note that although AKTW did not impose all of the 
many possible experimental constraints on the MSSM pa- 
rameter space, the resulting large model sample should 
be representative of the range of somewhat more difficult 
possibilities that one might face at the ILC. 

On our way to addressing the Inverse Problem at the 
ILC, we face the more immediate, and perhaps even more 
important, issue of the visibility of the various SUSY par- 
ticles in the AKTW models. We find that this is surpris- 
ingly non-trivial and is perhaps a more important task as 
one cannot differentiate between models which have no 
visible SUSY signatures. In our analysis below, we will 
perform a detailed study of the visibility of the various 
SUSY particles in all of the models. We will employ an 
extensive menu of search techniques and examine when 
they succeed and how they fail. Our philosophy will be to 
apply a general search strategy that performs uniformly 
well over the full MSSM parameter region, rather than 
make use of targeted searches for particular parameter 
points. We believe this mirrors the reality of an experi- 
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mental search for new physics and reflects the fact that 
not aU of the SUSY particles in these models will have 
been observed at the LHC (recall that the models we have 
inherited from AKTW are difficult cases at the LHC). ft 
is important to point out that given the large set of mod- 
els we examine, not linked to a particular SUSY breaking 
mechanism, provides a wide window on what studying 
generic SUSY may be like at the ILC. This is the first 
analysis of this kind and the results of this study were 
found to be unexpected and quite surprising, at least to 
the present authors. 

The possibility of measuring specific SUSY particle 
properties at the ILC for particular special points in 
the MSSM parameter space has a long history 0. Our 
approach here, however, provides several aspects which 
have not been simultaneously featured in earlier analy- 
ses: (i) We examine several hundred, essentially random, 
points in MSSM parameter space, providing a far wider 
than usual sampling of models to explore and compare. 
This gives a much better indication of how an arbitrary 
MSSM parameter point behaves and what experimental 
techniques are necessary to adequately cover the full pa- 
rameter space, (ii) We include all effects arising from 
initial state radiation (ISR), i.e., bremsstrahlung, as well 
as the specific ILC beamstrahlung spectrum for the su- 
perconducting RF design, including finite beam energy 
spread corrections. The beam spectrum is generated by 
GuineaPig |9[ [l^. (Hi) We incorporate all 2 ^ 2, 2^4 
and 2^6 SM background processes, including those 
resulting from initial state photons {i.e., from the corre- 
sponding 77 and 76^ interactions). These are generated 
with full matrix elements via WHIZ ARD/0 'MEGA Q 
for arbitrary beam polarization configurations and are 
fragmented using PYTHIA. There are well over 1000 of 
these processes |10 . (iv) We include ILC detector effects 
by making use of the java-based SiD Q detector fast 
simulation package org.lcsim |jl^, All in all, we be- 
lieve that we have performed our analysis in as realistic 
a manner as possible. 

Given the present cost and timing controversies in the 
world community surrounding the ILC (or, generally, any 
e~^e~ collider), this is a critical time to ascertain in as 
realistic way as possible the capabilities of such machines 
to explore new physics that may be discovered at the 
LHC. Thus, as will be described below, we present a very 
detailed study of the MSSM at the ILC using the SiD 
detector, an important milestone on the road to achieving 
this goal. 

The outline of this paper is as follows: Section 2 con- 
tains a discussion of the various kinematical features of 
the AKTW models under consideration, while Section 3 
provides an overview of our analysis procedure as well as 
a general discussion of the SM backgrounds. In the next 
Sections we separately consider the individual SUSY par- 



ticle analyses for sleptons (Section 4), charginos (Section 
5) and neutralinos (Section 6). This leads to an over- 
all set of model observation and comparisons in Section 
7 where we discuss the ability of the ILC to distinguish 
the AKTW models and resolve the LHC Inverse Prob- 
lem. This is followed by a discussion of our results and 
conclusions. 



SPECTRUM AND KINEMATICAL FEATURES 
OF THE MODELS 



Before beginning our analysis, we first examine the 
kinematical traits and features of the SUSY models that 
AKTW found to be indistinguishable at the LHC.@ 
This consists of a set of 383 models (i.e., 383 points in 
a 15-dimensional MSSM parameter space; we hereafter 
refer to distinct points in the MSSM parameter space 
as models). In their study, AKTW compared models 
pairwise, so that these 383 models correspond to 283 
pairs of models which gave indistinguishable signatures 
at the LHC. In some cases, models were found to give 
degenerate signatures multiple times. While this may 
naively seem to be a relatively small number of insepara- 
ble models, one needs to recall that AKTW performed a 
small sampling of a large parameter space (due to com- 
putational limitations). Based on the number of mod- 
els AKTW generated, the number of degeneracies they 
found led AKTW to estimate that a more complete sta- 
tistical sampling of the available parameter space volume 
would yield a degeneracy of each model with 0(10— 100) 
other points. 

One may wonder if there are any common features of 
these models that give rise to their indistinguishability at 
the LHC. AKTW demonstrated that these degeneracies 
are essentially the result of three possible characteristics 
that involve the relative composition of the physical elec- 
troweak gaugino sector in terms of the higgsino, wino, 
and bino weak eigenstates. These mechanisms are re- 
ferred to as 'Flippers', 'Sliders' and 'Squeezers' and are 
schematically shown in Figs. 0, and |[ The ambigui- 
ties that arise from these model characteristics originate 
directly from the manner in which SUSY is produced and 
observed at the LHC; several of these mechanisms can be 
simultaneously present. As is well-known, the (by far) 
dominant production mechanism for R-parity conserv- 
ing SUSY at the LHC is via the strong interactions, i.e., 
the production of squarks and gluinos. These particles 
then decay through a long cascade chain via the gener- 
ally lighter electroweak gaugino/higgsino partner states. 
This eventually leaves only the SM fields in the final state 
together with the stable Lightest SUSY Particle (LSP), 
which is commonly the lightest neutralino, appearing as 
additional missing energy. The decays of the SM fields 
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FIG. 1: Illustration of the Flipper ambiguity in the MSSM spectrum. The left panel displays a typical mass value for the 
Bino, Wino, and Higgsino mixing parameter in the weak eigenstate basis for two models A and B. The right panel shows the 
corresponding sparticle spectrum in two of the AKTW models, with the red arrows indicating the Flipper effect. 




FIG. 2: Illustration of the Slider ambiguity in the MSSM spectrum. The left panel displays a typical mass value for the 
Bino, Wino, and Higgsino mixing parameter in the weak eigenstate basis for two models A and B. The right panel shows the 
corresponding sparticle spectrum in two of the AKTW models, with the red arrows indicating the Slider effect. 
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FIG. 3: Illustration of the Squeezer ambiguity in the MSSM spectrum. The left panel displays a typical mass value for the 
Bino, Wino, and Higgsino mixing parameter in the weak eigenstate basis for two models A and B. The right panel shows the 
corresponding sparticle spectrum in two of the AKTW models, with the red arrows indicating the Squeezer effect. 



produce additional jets, leptons, and missing transverse 
energy from neutrinos. Unfortunately sleptons do not al- 
ways play a major role in these cascades (due to phase 
space considerations in the sparticle spectrum, see, e.g., 
the models in Fig. ^ so that much valuable information 
associated with their properties is generally lost. When 
comparing the possible decay chains which result from 
the produced squarks and gluinos, similar final states 
can occur if the identities of the higgsino, wino and bino 
weak states in the spectrum are interchanged while their 
masses are held approximately fixed. This is an exam- 
ple of the so-called 'Flipper' ambiguity (Fig. |^) where two 
spectra with interchanged electroweak quantum numbers 
can produce very similar final state signatures. A sec- 
ond possible source of degeneracy can arise from the fact 
that absolute masses, and in particular the mass of the 
LSP, are not well measured at the LHC in contrast to 
the mass differences between states |l5| . Thus models 
can have similar spectra but be somewhat off-set from 
each other in their absolute mass scale and hence be dif- 
ficult to distinguish; this represents the 'Shder' degener- 
acy (Fig. 1^). Lastly, pairs of states in the spectra with 
relatively small mass differences compared to the overall 
SUSY scale lead to relatively soft decay products in the 
cascade chain. Such a possibility can cause significant 
loss of information as well as general confusion in param- 
eter extractions and are termed 'Squeezers' (Fig. ||). Of 



course in all these cases some shifts are needed in the 
strongly interacting part of the SUSY spectrum to keep 
the various production rates and decay distributions com- 
parable between potentially indistinguishable models. It 
goes without saying that some degeneracies can also arise 
when more than one of these mechanisms are active si- 
multaneously. 

We now examine the physical particle spectra in the 
383 models found by AKTW to be indistinguishable at 
the LHC. First, we note that since AKTW have re- 
quired squarks and gluinos to have Lagrangian masses 
greater than 600 GeV in their parameter scans, the only 
states potentially accessible to the ILC will be the slep- 
tons and the sparticles associated with the electroweak 
gaugino/higgsino sector. Of particular phcnomcnologi- 
cal interest is the mass splitting between the Next-to- 
LSP(NLSP) and LSP (see Fig. |). Here, this is usually 
that between the lightest chargino, Xi , and lightest neu- 
tralino state, Xi- Generally this distribution for our set 
of AKTW models appears rather flat except for a huge 
and puzzling feature near ~ 270 MeV. It would seem 
that almost 40%, i.e., 141 of these models, experience 
this exact mass splitting between these two states. 

An investigation shows that this result is an artifact of 
the manner in which PYTHIA6.324 generates the phys- 
ical SUSY particle spectrum at tree-level from the La- 
grangian parameters. Recall that AKTW randomly gen- 
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TABLE I: Number of models at ^ = 500 GeV and 1 TeV 
which have a given final state kinematically accessible. Note 
that for the 500 GeV case, 96/242 models have only LSP or 
neutral pairs accessible while 61/242 models have no SUSY 
particles accessible. 



erated points in a 15-dimensional weak scale MSSM pa- 
rameter space, described in the previous Section, from 
which the physical SUSY particle masses are then cal- 
culated at tree-level via PYTHIA6.324. With this pro- 
cedure, it is possible that sometimes the mass of the 
lightest chargino xt turns out to be less than that of 
the Xi once the mass eigenstates are computed; this is 
usually considered to be 'unphysical' as it would imply 
charged Dark Matter in the standard cosmological pic- 
ture. PYTHIA6.324 handles this situation by artificially 
resetting the chargino mass to be greater than that of the 
LSP by Ti-i = TO^[)-|-2m7r without an associated warning 
message. This apparently happens frequently and causes 
the large peak in the distribution shown in Fig. ^. This 
feature is mentioned in the PYTHIA manual (where it is 
noted that the tree-level SUSY spectrum calculator is not 
for publication quality) , and has been further clarified in 

However, here we need 



later versions of PYTHIA |16 



to follow the analysis of AKTW as closely as possible 
to reproduce their sparticle spectra and specific model 
characteristics. Due to this and additional reasons dis- 



cussed below in the text, in our analysis we use a slightly 
modified version of PYTHIA6.324. In the strictest sense, 
these models are only 'unphysical' at the tree-level since 
loop corrections restore the correct mass hierarchy. We 
have checked that all 383 of the AKTW models have an 
appropriate mass spectrum when the SuSpect2.34 rou- 
tine ]l7[ |, which includes the higher order corrections, is 
employed to generate the physical spectrum. However, 
in the present work, the 141 models cannot be artificially 
saved simply by employing this mass re-assignment or 
by using SuSpect as their collider production and signa- 
ture properties would be modified as compared to the 
AKTW study. We thus drop them completely from fur- 
ther consideration. This leaves us with a sample of 242 
models which consist of 162 degenerate model pairs to 
examine. |Q 

Given these 242 models, we next address the question 
of what fraction of their SUSY spectra are kinematically 
accessible at a 500 or 1000 GeV ILC. The results are 
shown in Figs. ^ ^ and |^, which display the individ- 
ual mass spectra for the weakly interacting sectors of the 
various SUSY models under consideration. The full ac- 
cessible sparticle count for = 500 and 1000 GeV is 
presented in Fig. |^. There are many things to observe 
by examining these Figures. First, we recall from the 
discussion above that in all cases the squarks are too 
massive to be pair produced at the ILC so that we are 
restricted to the slepton and electroweak gaugino sec- 
tors. Here in Table | and in the Figures we see that for 
a 500(1000) GeV collider, there are only 22(137)/242, 
i.e., 22(137) out of 242, models with kinematically ac- 
cessible (which here means via pair production) selec- 
trons and smuons at ^/s = 500(1000) GeV; note that 
these two sparticles are degenerate in the MSSM. Sim- 
ilarly, 28(145)/242 of the models have accessible light 
staus, 6(55) of which also have kinematically accessi- 
ble selectrons/smuons. 53(92)/242 models have kine- 
matically accessible light charginos, 4(12) of which also 
have accessible selectrons/smuons and 6(12) of which also 
have accessible staus. At y/s = 1 TeV, 19 of these 92 
models with accessible light charginos also have the sec- 
ond chargino accessible by pair production. Very im- 
portantly, at ^/s = 500 GeV, in 96/242 models the only 
kinematically accessible sparticles are neutral, e.g., Xi 
or i>, while 61/242 other models have no SUSY particles 
accessible whatsoever. At y/s = 1 TeV, these numbers 
drop to only 0/242 and 3/242 in each of these latter cat- 
egories, respectively. Recalling that wc are looking at es- 
sentially random points in the MSSM parameter space, 
we see from this simple counting exercise that 60% 
of the models will have no 'traditional' SUSY signatures 
at a 500 GeV ILC, whereas a 1 TeV machine essentially 
covers almost all the cases. This is a strong argument 
for having the capability of upgrading to 1 TeV at the 
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FIG. 4: Low mass end of the lightest chargino-LSP mass difference which displays the PYTHIA feature in the sparticle spectrum 
generator discussed in the text. 



ILC as quickly as possible. However, in the analysis that 
follows we will consider only the case of a 500 GeV ILC 
with the 1 TeV case to be considered separately in the 
future. Table ^ summarizes the kinematic accessibility of 
all the relevant MSSM final states for y/s = 500 GeV in 
our study as well as the corresponding results for 1 TeV. 

Given that so many models have such a sparse SUSY 
spectrum at = 500 GeV, it is not uncommon for one 
of the two models in the pair we are comparing to have 
no kinematically accessible sparticles. In such a case, 
breaking the model degeneracy at the LHC might seem 
to be rather straightforward, as for one model we might 
observe SUSY signals above the SM background but not 
for the other in the pair. Of course, at the other end of 
the spectrum of difficulty, one can imagine cases where 
both models being compared are Squeezers, in which case 
model differentiation will be far more difficult and having 
an excellent ILC detector will play a much more impor- 
tant role. 
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FIG. 5: Selectron/smuon (top) and stau (bottom) mass spectra in the 242 models under study. 
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FIG. 6: Sneutrino mass spectra in the 242 models under study. 
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FIG. 7: Chargino (top) and neutraiino (bottom) mass spectra in the 242 models under study. 
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FIG. 8: Number of models with various kinematically accessible MSSM states at a 500 (red) or 1000 (blue) GeV e^e collider 
via pair production. 
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ANALYSIS PROCEDURE AND GENERAL 
DISCUSSION OF BACKGROUND 

To determine whether or not the ILC resolves the LHC 
inverse problem, we compare the ILC experimental sig- 
natures for the pairs of SUSY models that AKTW found 
to be degenerate, and see whether these signatures can be 
distinguished. We examine numerous production chan- 
nels and signatures for supersymmetric particle produc- 
tion in e^e~ collisions. Before the model comparisons 
can be carried out, we must first ascertain if the pro- 
duction of the kinematically accessible SUSY particles is 
visible above the SM background. Our analysis proce- 
dure is described in this Section. 



Event Generation of Signal and Background 

We generate 250 fb^i of SUSY events at ^/s = 
500 GeV for each of the AKTW models for both 
80% left- and 80% right-handed electron beam po- 
larization with unpolarized positron beams, pro- 
viding a total of 500 fb~^ of integrated lumi- 
nosity. To generate the signal events, we use 
PYTHIA6.324 g in order to retain consistency with 
the AKTW analysis. However, as will be described in 
detail below, we find that PYTHIA underestimates the 
production cross section in two of our analysis channels, 
and in these two cases we employ CompHEP We 
also analyze two statistically independent 250 fb^^ sets of 
Standard Model background events for each of the two 
electron beam polarizations. We then study numerous 
different analysis channels. When wc determine if a sig- 
nal is observable over the SM background in a particular 
channel, we statistically compare the combined distribu- 
tion for the signal plus the background from our first 
background sample with the distribution from our sec- 
ond, independent background sample. When we perform 
the model comparisons, we add each set of SUSY events 
to a distinct Standard Model background sample gener- 
ated for the same beam polarization. We then compare 
observables for the many different analysis channels for 
these two samples of signal and background (i.e., model 
A + background sample 1 is compared to model B -|- 
background sample 2). It is important to note that we 
take into account the full Standard Model background in 
all analysis channels rather than only considering the pro- 
cesses that are thought to be the dominant background to 
a particular channel; surprisingly, sometimes many small 
contributions can add up to a significant background. 

Our background contains all SM 2 ^ 2, 2 ^ 4, and 
2 — !■ 6 processes with the initial states e^e~, e^7, or 
77; in total there are 1016 different background chan- 
nels. These events were generated by T. Barklow 110] 



with O'MEGA as implemented in WHIZARD |l|, which 
uses full tree-level matrix elements and incorporates a 
realistic beam treatment via the program GuineaPig [|j . 
The use of full matrix elements leads to qualitatively dif- 
ferent background characteristics in terms of both to- 
tal cross section and kinematic distributions compared 
to those from a simulation that uses only the produc- 
tion and decay of on-shell resonances, e.g., the procedure 
generally employed in PYTHIA. WHIZARD models the 
flux of photons in e^7 and 77 initiated processes via the 
equivalent photon approximation. However, in the stan- 
dard code, the electrons and positrons which emit the 
photon(s) that undergo hard scattering do not receive a 
corresponding kick in px, in contrast to the electrons or 
positrons that undergo initial state radiation. The ver- 
sion of WHIZARD used here to generate the background 
events was thus amended to correct this slight inconsis- 
tency in the treatment of transverse momenta. An il- 
lustration of the resulting effects from employing exact 
matrix elements and modelling the transverse momentum 
distributions in a realistic fashion is presented in Fig. ||. 
This Figure compares the transverse momentum distri- 
bution for the process e^e~ e^e~i'ei^e in the SM after 
our selectron selection cuts (see Section 4.1) have been 
applied, as generated with PYTHIA versus the modified 
version of WHIZARD, using the same beam spectrum in 
both codes. We see that in this case, the px distribution 
generated by PYTHIA is smaller and has a shorter tail. 

We now discuss our treatment of the beam spectrum 
in further detail. The backgrounds were generated us- 
ing a realistic beam treatment, employing the program 
GuineaPig 1^ to model beam-beam interactions. Finite 
beam energy spread was taken into account and com- 
bined with a beamstrahlung spectrum specific to a cold 
technology linear collider, i.e., the ILC. The effect of 
beamstrahlung is displayed in Figure |l^, which shows 
the invariant mass of muon pairs formed by e+e" colli- 
sions with the beam spectrum we employ. The result- 
ing spectrum is somewhat different qualitatively from 
a commonly used purely analytic approximate approach 

While our backgrounds contain 500 fb^^ of total inte- 
grated luminosity for processes with initial e^e~ or e^7 
states, some 77 initiated processes yield very high cross 
sections, and thus a smaller number of events had to be 
generated and then rcscalcd due to limited storage space. 
In total, our background sample uses approximately 1.7 
TB of disk space. This rescaling of some 77 processes in- 
troduces artificially large fluctuations in the correspond- 
ing analysis distributions. In order to remedy this, we 
employ the following procedure: we combine the two in- 
dependent background samples for the affected reactions, 
and then randomly reallocate each entry on a bin by bin 
basis to one of the two background sets. Thus, on aver- 
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visible pT, selectron analysis, e- = -80% pol. 
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■ total PYTHIA (WW, gamma) 
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FIG. 9: Transverse momentum distribution in e+e^ e^e~VgUf. as generated via PYTHIA and WHIZARD, for 250 fb~^ 
of integrated luminosity with 80% left-handed electron beam polarization at ^/s = 500 GeV. Our selectron selection cuts 
(discussed in the text below) have been applied. 
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FIG. 10: Invariant mass of the muon pairs in e^e — > ^'^ ^ at -/s = 500 GeV, using the beamspectrum described in the text. 
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age, each histogram contains an equal amount of entries 
bin per bin, while remaining statistically independent. Of 
course, this procedure does not completely eliminate the 
fluctuations. However, due to the random reallocation of 
entries, the contribution of these fluctuations to the sta- 
tistical analyses in our comparison of models performed 
in Section 7 is greatly reduced. 



Analysis Procedure 

For each SUSY production process, we perform a cuts- 
based analysis, and histogram the distributions of various 
kinematic observables that we will describe in detail in 
the following Sections. We apply a general analysis strat- 
egy that performs uniformly well over the full MSSM pa- 
rameter region. Each analysis is thus applied to every 
model in exactly the same fashion; there are no free pa- 
rameters, and we do not make use of any potential infor- 
mation from the LHC; in particular we assume that the 
LSP mass is not known. Recall that the AKTW models 
that we have inherited are difficult cases at the LHC, and 
thus in general we cannot make any assumptions about 
what measurements, if any, will have been performed by 
the LHC detectors. We also note that AKTW did not 
impose any additional constraints from flavor physics, 
cosmological observations, etc. Such a global analysis 
is clearly desirable but is beyond the scope of the present 
study and is postponed to a future publication. 

Our background and signal events described above 
are piped through a fast detector simulation using the 
org.lcsim detector analysis package [^], which is cur- 
rently speciflc to the SiD detector design p2[. o rg.lcsim 
is part of the Java Analysis Studio (jas3)[E2|, a gen- 
eral purpose java-based data analysis tool. The org.lcsim 
fast detector simulation incorporates the specific SiD de- 
tector geometry, finite energy resolution, acceptances, as 
well as other detector specific processes and effects. Un- 
fortunately, the identification of displaced vertices and a 
measurement of dE/dx arc not yet implemented in the 
standard, fully tested version of the simulation package 
employed here, although preliminary versions of these 
functions are under development. The present study rep- 
resents the first large end user application of the Icsim 
software package and hence we prefer to use the standard, 
tested, version of the software without these additional 
features. The output of our Icsim-based analysis code is 
given in terms of AID A histograms, where AID A refers 
to Abstract Interfaces for Data Analysis |^ , a standard 
set of Java and C-f + interfaces for creating and manip- 
ulating histograms, which is incorporated into the jas3 
framework. 

org.lcsim allows for the study of various different de- 
tectors, whereby an xml description of the specific de- 



tector is loaded into the software in a modular fashion. 
Currently, xml descriptions for various slightly differ- 
ent versions of the SiD detector geometry are publicly 
available. We use the SiD detector version studied ex- 
tensively at Snowmass 2005 (sidaugOS) [Q. In addi- 
tion, two files, called ClusterParameters .properties 
and TrackingParameters .properties allow the user to 
adjust various tracking and energy resolution parameters. 
We set the parameters such that we closely follow the SiD 
detector outline document (DOD) jl^. In particular, we 
employ the following configuration in our study: 

• The minimum transverse momentum of registered 
tracks is given by pt > 0.2 GeV. 

• There is no tracking capability below 142 mrad, 
which corresponds to \cos9\ < 0.99. 

• Between 142 mrad and 5 mrad, electromagnetically 
charged particles appear as neutral clusters. 

• There is no detector coverage below 5 mrad. 

• The jet energy resolution is set to 30% /Ve. 

• The electromagnetic energy resolution is set to 
18%/Ve. 

• The hadronic energy resolution is set to 50% /Ve. 

• The hadronic degradation fraction is r = 1.0. 

• The electromagnetic jet energy fraction is Wj = 
0.28. 

• The hadronic jet energy fraction is Wh = 0.1. 

For a detailed explanation of these parameters we refer 
to the SiD DOD, specifically Section IV. B regarding the 
energy resolution parameters. 

However, we note that the lack of tracking capability 
below 142 mrad causes highly energetic forward muons 
to not be reconstructed. They are too energetic to de- 
posit energy into clusters and are thus undetected and 
appear as missing energy. This effect produces a sub- 
stantial Standard Model background to, e.g., our stau 
analysis (see Section ), where we allow one tau to de- 
cay hadronically and the other leptonically. In this case, 
we keep events with one electron and one muon of op- 
posite charge, which can be mimicked by 77 — > 
events where one of the beam electrons is kicked out 
sufficiently to be detected, but one of the final state 
muons is too energetic and too close to the beam axis 
to be reconstructed. We find that this background is 
substantial and, given these detector parameters, can 
only be eliminated by discarding all tau events with elec- 
trons/positrons in the final state. 
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The default jet finding algorithm of org.lcsim is the 
JADE jet algorithm in the E scheme with y^^t = 
0.005 employed as the default setting. The JADE jet 
algorithm in the E scheme is defined as follows: 

mm{pi +Pj)'^ = mm2EiEj{l - cos%) > y^^tS (1) 
Pij = Pi + Pj for the recombination scheme 

The default y^^, however, is too small, and causes soft 
gluons to produce far too many jets. We therefore set the 
value of t/cut to t/^ut = 0.05 within the JADE jet algorithm. 
In addition, one must take care when using the default 
org.lcsim jet finder, as every parton, including leptons 
and photons, is in principle identified as a jet. More 
sophisticated jet finders are in the development stage. 
We thus use the default jet finder, but with additional 
checks on the jet particle content to discard non-hadronic 
"jets". 

We perform searches for slepton, chargino, and neu- 
tralino production and in many cases design analyses for 
several different decay channels of these sparticles. Each 
of our analyses is designed to optimize a particular signa- 
ture, and we apply each analysis to every AKTW model. 
A particular model may or may not produce a visible sig- 
nature in a specific channel. We will describe our cuts in 
detail for each analysis channel in the Sections below. 

As a starting point, we incorporate sets of kinematic 
cuts that were developed in various previous supersym- 
metric studies in the literature (the specific references are 
given in the following Sections). However, in many cases 
we find that some of these cuts are optimized for spe- 
cific Supersymmetry benchmark points, e.g., the Snow- 
mass Points and Slopes (SPS) and are too stringent 
for the general class of models we study here. In other 
cases, we find that the cuts employed in the literature 
are not stringent enough to sufficiently reduce the Stan- 
dard Model background in order to obtain a good sig- 
nal to background ratio. Through a seemingly endless 
series of iterations, we have thus designed sets of cuts 
(which are described in the following Sections) that opti- 
mize the signal to background ratio for an arbitrary point 
in MSSM parameter space. This scenario corresponds to 
a first sweep of ILC data in search of a SUSY signal, 
and is therefore a reasonable course of action. We also 
remind the reader that these AKTW models are difficult 
at the LHC and hence the slepton, chargino, and neu- 
tralino states will not necessarily be observable at the 
LHC. 

The signatures that we have developed analyses for are 
summarized in Table which lists the signature, dom- 
inant background source, and the observable kinematic 
distribution for each SUSY production process. We note 
that in some cases, the same signature can arise from dif- 
ferent sources of sparticle production, e.g., fi'^fi~+ miss- 
ing energy can occur from both smuon and chargino pro- 



duction. Indeed, it is well known that sometimes SUSY 
is its own background and we will note this in the fol- 
lowing Sections. Our cuts, however, are chosen such to 
minimize this effect. 

As discussed in the introduction, the first step in our 
analysis is to determine whether or not a given SUSY 
particle is visible above the SM background. Specifically, 
for a kinematic distribution resulting from our analysis 
of a given observable, we ask whether or not there is suf- 
ficient evidence to claim a 'discovery' for a SUSY particle 
within a particular model. There are many ways to do 
this, but we follow the Likelihood Ratio method, which 
we base on Poisson statistics. (See, e.g., Ref. |2^). In 
this method, we introduce the general Likelihood distri- 
bution: 



bins 



(2) 



where ni{fj,i) are the number of observed (expected) 
events in each bin i and we take the product over all 
the relevant bins in the histogram. As discussed above, 
we have generated two complete and statistically inde- 
pendent background samples, which we will refer to as 
Bl and B2. Combining the pure signal events, S, which 
we generate for any given model with one of these back- 
grounds, we form the Likelihood Ratio 



B = L{S + Bl, B2)/L{B1, B2) . 



(3) 



The criterion for a signal to be observed above back- 
ground is that the significance, S, satisfy 



S = V21og^ > 5. 



(4) 



This corresponds to the one-sided Gaussian probability 
that a fluctuation in the background mimics a signal of 
~ 2.9 • 10~^, which is the usual ha discovery criterion. 
When employing this method, we sometimes encounter 
bins within a given histogram for which there is no back- 
ground due to low statistics but where a signal is ob- 
served. In this case, the function L is not well-defined. 
When this occurs we enter a single event into the empty 
background histogram in that bin. 

Given that our full SM background samples are only 
available at fixed energies, our toolbox does not include 
the ability to do threshold scans. As is well known, this 
is a very powerful technique that can be used to ob- 
tain precision mass determinations for charged SUSY (or 
any other new) particles that are kinematically accessi- 
ble. Such measurements would certainly aid in the dis- 
crimination between models, especially in difficult cases 
where the measurements we employ do not suffice. In 
addition, especially for sparticles which decay inside the 
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Sparticles Produced 


Signature 


Main Background 


Observable 




e^e~ + missing E 




energy of e~ 




ii^ ii~ + missing E 


e+e" ji^vivi^j^ , 


energy of /i^, /i~ 


T+T- 


t^t" + missing E 


e*7 —> e^vu 


energy of tau jets 






e±7 ^ e^l+l- 




UV* 


jjjjl^l' + missing E 


77 — > cc, hb 


missing energy 


vv* 


+ nussmg E 


none 


missing energy 




+ missing E 




energy of , ji" 


xXxi 


jjfj,^ + missing E 


e^e~ qq'li'i 


energy and invariant 






77 — > T^T~ 


mass of dijet pair 






77 qq 




xXxl 


jjjj + missing E 


77 qq 


energy and invariant 








mass of dijet pairs, 








missing energy 


xXxi 


7 + charged tracks 




recoil mass 






e±7 ^ e^l+r 




xXxl or 


2 stable charged tracks 


e+e- e+e" +ISR,BS 


p/E 




l^l^ + missing E 


77 -> i+r 


invariant mass 








of lepton pair 


X^X? 


jj + missing E 


e^-y Veqq 


invariant mass 








of jet-pair 


x?x? 


7 + nothing 




photon energy 



TABLE II: Summary of signatures and observables in all analysis channels that we study and sources of the main standard 
model background. I = e, 



detector volume, input from the excellent SiD vertex de- 
tector could prove extremely useful. In the analysis pre- 
sented below, the vertex detector is used only for track 
matching and not as a search tool for long-lived states. 



SLEPTON PRODUCTION 
Charged Slepton Pair Production 

For detecting the production of charged sleptons, we 
focus on the decay channel 

l+t ^l+l-X°iX°i, (5) 

that is, the signature is a lepton pair plus missing energy. 
In the cases of selectrons and smuons these signatures are 
fairly straightforward to study; the stau case is slightly 
more complicated due to the more involved tau identifi- 
cation. 

As is well known, the main Standard Model back- 
ground for all of these cases arises from the produc- 
tion of W pairs followed by their subsequent decay into 
lepton-neutrino pairs and from Z-boson pair production, 
where one Z decays into a charged lepton pair and the 
other into a neutrino pair. A significant background also 



arises from gamma-induced processes through beam- and 
bremsstrahlung. 

The W pair background produces leptons that are pre- 
dominantly along the beam axis towards |cos0| « 1, 
where 9 takes on the conventional definition. This is 
because the decaying W bosons are produced either 
through s-channel Z- or 7-exchange, for which the dif- 
ferential cross section is proportional to (1 -I- cos^O), or 
through t-channel neutrino-exchange, which behaves as 
1/ sin^(6'/2). The photon-induced background also yields 
electrons that are peaked along the beam axis because 
they are mainly produced at low px from beam- and 
bremsstrahlung, although our more realistic beamspec- 
trum has a larger pr tail than the PYTHIA-generated 
backgrounds studied conventionally (cf. the discussion in 
Sec. ). As we will illustrate below in Section 4.1.3, hav- 
ing the best possible forward detector coverage in terms 
of tracking and particle identification (ID) is therefore of 
utmost importance to reduce the Standard Model back- 
ground. 

To reduce the SM background, we employ a series of 
cuts that have been adapted and expanded from previous 
studies ||2^, ^ |2^. Our cuts are fairly similar for all 
slepton analyses. We will discuss them in detail in our 
selectron analysis presented below, and then will list the 
cuts with only brief comments in our discussion of smuon 
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and stau production. 

Selectrons 

As discussed above, in the case of selection production 
we study the clean decay channel 

e+e- ^ e+e-x?X? , (6) 

that is, the signature is an electron pair plus missing 
energy. The main backgrounds arising from the SM orig- 
inate in W and Z pair production, followed by their 
leptonic decays, along with several processes originat- 
ing from both 77 and 76 interactions. To reduce these 
backgrounds we employ the following cuts, which are ex- 
panded from those in p7| : 

1. We require exactly two leptons, identified as an 
electron and a positron, in the event and that there 
be no other charged particles. This removes SM 
backgrounds where, for example, both Z bosons 
decay into charged leptons. 

2. E^,, < 1 GeV for | cos 6*1 > 0.9, where cor- 
responds to the visible energy in the event. This 
helps to reduce large SM background from forward 
W production, as well as beam-/bremsstrahlung re- 
actions that yield leptons predominantly along the 
direction of the beam axis. 

3. < OAy/s in the forward hemisphere. Here, the 
forward hemisphere is defined as the hemisphere 
around the thrust axis which has the greatest vis- 
ible energy. Since we only have 2 visible parti- 
cles in the final state, this amounts to defining 
the forward hemisphere about the particle with the 
highest energy. The SUSY signal has missing en- 
ergy in both hemispheres, whereas the SM reaction 
e^e~ — *■ ZZ — > e'^e~vv has missing energy in only 
one of the hemispheres since the decay Z ^ vv oc- 
curs in the hemisphere opposite of the Z decay to 
charged leptons. 

4. The angle between the reconstructed electron- 
positron pair is restricted to have cos0 > —0.96. 
Since SUSY has a large amount of missing energy, 
the selection pair will not be back-to-back, in con- 
trast to the SM background events. 

5. We demand that the visible transverse momentum, 
or equivalently, the transverse momentum of the 
electron-positron pair, pTvis — Vt ^ 0.04y^. 
This cut significantly reduces both the 77 and 6^7 
backgrounds which are mostly at low px- 



6. The acoplanarity angle must satisfy Acj)'^ ^ > 40 
degrees. Since we demand only an electron and 
positron pair, the acoplanarity angle is equivalent 
to TT minus the angle between the transverse mo- 
mentum of the electron and positron, Acjj'^ — 
n — 9t- This requirement translates to a restriction 
on the transverse angle of cos9t > 0.94. This cut 
further reduces contributions from both the W^-pair 
and 77 backgrounds where the e+e^ pair tends to 
be more back-to-back. 

7. A4+e- < Afz - 5 GeV or M^+^- > Mz + 5 GeV, 
where Afg+g- is the invariant mass of the lepton 
pair. This cut is to further remove events from Z 
boson pair production with subsequent decays into 
e^e~ pairs. 

As already mentioned above, we note that below 142 
mrad (|cos0| > 0.99), the SiD detector does not have 
particle tracking information according to the current de- 
tector design , and any charged particle in this region 
appears only as a neutral electromagnetic cluster. How- 
ever the first cut listed above, where we demand exactly 
an e^e~ pair in the final state, substitutes for poten- 
tially more detailed cuts that assume tracking capabili- 
ties down to much smaller angles. 

The standard selectron search analysis is based on the 
energy distribution of the final state electron or positron. 
Since the selectron decays into a clean 2-body final state, 
the e~/e+ energy distribution has a box-shaped "shelf" 
in a high statistics, background- free, perfect detector en- 
vironment in the absence of radiative effects. Kinematics 
determines the minimum and maximum electron energies 
which are related to the two unknown masses of the se- 
lectron and LSP by 

The sharp edges of this box-shaped energy distribution 
allow for a precise determination of the selectron and 
LSP masses. However, due to beamstrahlung, the ef- 
fective above will vary, and once detector effects are 
also included the edges of this distribution will tend to 
be slightly washed out. Since our goal here is to simply 
detect superpartners and then distinguish models with 
different sets of underlying parameters, we do not need 
to perform a precise mass determination in the present 
analysis. We also consider additional kinematic observ- 
ables, such as the distribution of p^'^ and the e+e^ in- 
variant mass Afgg, as they will be useful at separating 
different SUSY signal sources. 
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The successive effect of each of the above cuts on the 
SM background is iUustrated in Fig. ^ Here, we show 
the electron and positron energy distribution for 250 fb^^ 
of simulated Standard Model background for RH electron 
beam polarization at yfs = 500 GeV. The y-axis corre- 
sponds to the number of events per 2 GeV bin. We note 
that cuts number 1-5 essentially eliminate any potential 
background arising from the large Bhabha scattering and 
77 — + e+e^ cross sections. The main contribution to 
the background remaining after these cuts arises from 
processes involving W and Z pair production from elec- 
tron positron initial states, i.e., e+e^ t'lifViil'^ ^ with 

) = e, T, as shown in Fig. |l2|. Wc find that most of the 
photon initiated background has been removed by our 
cuts. Note that applying these cuts in a different order 
would necessarily show a different level of effectiveness. 

A further comment on cut number 5 is in order. One 
finds that increasing this cut from 'p^ ^ > 20 GeV to 
Vt^~ > 30 GeV at = 500 GeV to further reduce 
the photon initiated background, introduces a dip in the 
center of the "shelves" for MSSM models that have edges 
near E^, ~ 30-40 GeV. This apparently occurs when both 
visible leptons have approximately the same amount of 
energy, so that their transverse momenta partially cancel, 
leaving insufficient visible transverse momentum to pass 
the cut. When one of the leptons is more energetic than 
the other, the visible pt is generally above the cut. Thus 
increasing the cut on visible pr in order to reduce the 
background further also affects the signal in a perhaps 
somewhat unexpected way. The same observation also 
applies in the smuon analysis below. 

Figure |l^ shows how a typical signal from a model 
with kinematically accessible selectrons responds to the 
same cuts that were applied to the backgrounds above. 
Note that while the cuts reduce the backgrounds by many 
orders of magnitude, the signal is reduced only by a factor 
of 2 - 3. 

We now examine selectron pair production for the 
AKTW models. In this case, there are 22/242 (22 out 
of 242) models with kinematically accessible selectrons 
at ^/s = 500 GeV. The 6^(7^) is accessible in 9(15) of 
these 22 models; for 2 models we find that both states 
are potentially visible. Note that fewer than 10% of our 
models have this relatively clean channel accessible. Se- 
lectrons are pair produced via s-channel 7 and Z ex- 
change as well as i-channcl neutralino exchange. For the 
case of the well-studied e ex? decay mode, which we 
examine here, selectrons are usually searched for by ex- 
amining the detailed structure of the resulting individual 
e^ energy spectra and looking for any excesses above 
the expected SM background. As is by now well-known 
and briefly discussed above, in the absence of such back- 
grounds, with high statistics and neglecting any radiative 
effects, the 2-body decay of the selectrons lead to flat, 



horizontal shelves in this distribution. In a more realistic 
situation where all such effects are included and only fi- 
nite statistics are available, the general form of the shelf 
structure remains but they are now jagged, tilted down- 
wards (towards higher e^ energies), and have somewhat 
smeared edges. These effects are illustrated in Fig. |l^ 
which shows examples of the e^ spectra (adding signal 
and background) for some representative AKTW models 
containing kinematically accessible selectrons with either 
beam polarization configuration. There are several im- 
portant features to note in these Figures. The detailed 
nature of the e signal in the energy spectrum shows 
significant variation over a wide range of both magni- 
tude and width depending upon the Sl.r and Xi masses 
and the resulting production cross sections. Recall that 
i-channel sneutrino exchange can be important here and 
dramatically affects the size of this cross section. The 
ratio of signal to background is not always as large as 
in most cases discussed in the literature. In addition, 
we note that the range of possible signal shapes relative 
to the SM background can be varied; not all of our sig- 
nals appear to be truly shelf-like. In some models, the 
background overwhelms the signal. Note that RH polar- 
ization leads to far smaller backgrounds than does LH 
polarization as would be expected, this being due to the 
diminished contribution from VF-pairs which prefer LH 
coupling. 

Of the 22 kinematically accessible models, 18(15) lead 
to signals with a visibility significance over background of 
5 > 5 at these integrated luminosities assuming RH(LH) 
beam polarization. Combining the LH and RH polar- 
izations channels we find that 18/22 models with selec- 
trons lead to signals with significance > 5. Furthermore, 
8/9 models with kinematically accessible cl are observ- 
able while 12/15 models with en are visible. Note that 
4 models have selectrons with masses that are in excess 
of 241 GeV. This leads to a strong kinematic suppres- 
sion in their cross sections and, hence, very small signal 
rates, so they are missed by the present analysis. Some 
of the models in both the RH and LH polarization chan- 
nels have a rather small S/B and are not easily visible 
at this level of integrated luminosity; typical examples of 
'difficult' models are presented in Figs. |l^ and 

It is interesting to compare these results to what we 
obtain in the case of the well-studied SPSla' benchmark 
model [^; Fig. [l^ shows the electron energy distribution 
for this model for both beam polarization choices in an 
analogous manner to that shown in the previous Figures 
for the AKTW models. Due to the large production cross 
section, detecting the signal in this case is rather trivial 
as we would expect from the detailed studies made in the 
literature. The most important thing to notice from this 
Figure is that SPSla' leads to substantially larger signal 
rates than in any of the models we are investigating in 
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the present analysis. In fact, the SPSla' signal rates can 
be almost two orders of magnitude larger than some of 
the models we are examining here. We also observe the 
obvious presence of two shelves, especially in the case of 
LH beam polarization, clearly indicating that both the 
and e_R states are kinematically accessible and are being 
simultaneously produced. 

Interestingly, one finds that there are a number of mod- 
els, particularly in the case of RH polarization, which do 
not have kinematically accessible selectrons but which 
have visible signatures in the e-pair analysis. This is an 
example of SUSY being a background to SUSY. There 
are, of course, other SUSY particles which can decay 
into and missing energy, e.g., chargino pair produc- 
tion followed by the decay xt W^*Xi with W* — > ev, 
or associated X2X1 production followed by the decay 
X2 e'''e~Xi- Both of these processes result in the same 
observable final state. Figure |l^ shows some of these 
'fake' models that appear in our selectron analysis in the 
case of RH polarization. Fake models, by which we mean 
models where other SUSY particle production leads to 
a visible signature in the selectron analysis, also appear 
for LH polarization and, in fact, we find 14 counterfeit 
models for either polarization. Note that the shapes of 
these fake model signatures are somewhat different than 
those in a typical model with actual selectrons present; 
there are no truly shelf-like structures and the ener- 
gies are all peaked at relatively low values. This occurs 
because in these examples the final state electrons are the 
result of a 3- (or more) body decay channel when the W 
boson is off-shell and because the xf ~ Xi mass splitting 
is relatively small. Both of these conditions are present 
in most of our models. 

If further differentiation from the fake signatures is re- 
quired, we need to examine a different kinematic distribu- 
tion, e.g., the invariant mass of the e+e^ pair, Mee. One 
would expect the counterfeit signals to populate small 
values of Alee while the models with actual selectrons will 
have a higher number of events with larger values of Mee- 
This is indeed the case as can be clearly seen in Fig. |l^. 
As we will sec below, fake signals occur quite commonly 
in almost all of our analyses. Though specific analyses 
are designed to search for a particular SUSY partner it is 
quite easy for other SUSY states to also contribute to a 
given final state and be observed instead, e.g., a similar 
signature can be generated using the visible px of the 
electron. 

It would be interesting to return to this issue with a 
wider set of models that lead to larger mass splittings in 
the electroweak gaugino sector to see how well selectrons 
and charginos can be differentiated under those circum- 
stances. 
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FIG. 11: Standard Model background after each of the cuts listed in the text is successively imposed for an incoming right- 
handed electron beam with 80% polarization and 250 fb~^ of integrated luminosity. 
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FIG. 12: Remaining SM background after the full set of selectron selection cuts listed in the text are imposed. This is 
generated from 250 fb~^ of SM events with 80% right-handed (solid blue line) and 80% left-handed (solid red line) electron 
beam polarization, and unpolarized positron beam at ^/s — 500 GeV. The dotted lines show the main process contributing to 
the background, e'^e^ — > £~ uiUiii'^ , with S = e,T, for 80% right-handed (dotted pink line) and 80% left-handed (dotted 
green line) electron polarization. Note that here and in other Figures below the spikes in the full background and the main 
contributions are because of rescaling issues and the thus necessary mixing of the two independent background samples (cf. 
the discussion in Section 3.1). This mixing is a random procedure which explains why the spikes are not all in the same bins. 
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FIG. 13: The electron energy distribution for selectron production after successively imposing each of the cuts listed in the 
text for the case of a right-handed incoming electron beam with 80% polarization and 250 fb~^ of integrated luminosity. 
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FIG. 14: Electron energy distribution: the number of events/2 GeV bin after imposing the full set of cuts discussed in the text 
for several representative models. RH(LH) beam polarization is employed in the top(bottom) panel, assuming 80% electron 
beam polarization an integrated luminosity of 250 fb~^ for either polarization. The SM background is shown as the black 
histogram. 
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FIG. 15: Same as the previous Figure but now for three models which are difficuh to observe due to small cross sections in 
RH polarization channel. 
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FIG. 16: Same as the previous Figure but now for the LH polarization channel for two sets of models which are a bit more 
difficult to observe due to small cross sections. 
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FIG. 17: The electron energy distribution after imposing the lull set ol cuts discussed in the text for the benchmark model 
SPSla'. RH(LH) beam polarization is employed in the top(bottom) panel, assuming an integrated luminosity ol 250 fb~^ for 
either polarization. The SM background is shown as the black histograms. 
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FIG. 18: The electron energy distribution after imposing the lull set ol cuts discussed in the text lor a subset ol lake selectron 
models in the RH polarization channel. The SM background is shown as the black histogram as usual. 
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FIG. 19: The dielectron invariant mass spectrum in the RH polarization channel for true selectron models (top) and those 
caused by other SUSY particles (bottom). 
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Smuons 

For jl pair production the standard search/ analysis 
channel is 

fi+r ^ A*+/i-x?X? , (9) 

i.e., the signature is a muon pair plus missing energy. 
Smuon production occurs via s-channel 7 and Z ex- 
change; there is no corresponding t-channel contribution 
as in the case of selectrons. As in the selection analysis, 
the dominant background arises from leptonic decays of 
VF-pair and Z-pair production, as well as the ubiquitous 
(though somewhat less important in this case) 77 back- 
ground. Since the background and signal are similar to 
those for selectron production, our choice of cuts here will 
follow those employed in the selectron analysis above and 
are adapted from those proposed by Martyn |^ (see also 

1. No electromagnetic energy (or clusters) > O.Oly^ 
in the region | cosO\ > 0.995. 

2. Exactly two muons are in the event with no other 
charged particles and they are weighted by their 
charge within the polar angle —0.9 < Q^cosO^ < 
0.75 with no other visible particles. This removes 
a substantial part of the ly-pair background. 

3. The acoplanarity angle satisfy A(/)^^ > 40 degrees. 
This reduces both the VK-pair and 77 backgrounds. 

4. |cos6lp„„_J < 0.9. 

5. The muon energy is constrained to be > 
0.004Vi. 

6. The transverse momentum of the dimuon system, 
or equivalently, visible transverse momentum (since 
only the muon pair is visible), satisfy pxvis = Pt^ > 
0.04-^5. This removes a significant portion of the 
remaining 77 and 6^7 backgrounds. 

The remaining SM background after these cuts have 
been imposed are displayed in Fig. |2^ for both polar- 
ization configurations. The main background to /i-pair 
production, e~^e~ l~uivirl'^ , I = fi^T, is also shown in 
the Figure, and we see that it essentially comprises the 
full background sample. The background is somewhat 
smaller here than in the case for selectron production, 
as beam remnants from 7-induced reactions are not con- 
fused with the signal for smuon production. 

As in the case of selectron production, there are only 
22 out of 242 models which have kinematically accessi- 
ble smuons at ^/s = 500 GeV. The /iL(_R) is found to 
be kinematically accessible in 9(15) of these cases, there 



again being 2 models where both smuon states can be si- 
multaneously produced. In the jl fiXi decay channel, 
smuons are observed by detecting a structure above the 
SM background in the muon energy distribution, similar 
to the search for selectrons. Here too, as is well-known, 
in the absence of such backgrounds, with high statistics 
and neglecting radiative effects, the 2-body decay of the 
/2's leads to horizontal shelf-like structures. In the more 
realistic situation where all such effects are included, the 
shelves remain but are now tilted downward (towards 
higher muon energies), as in the selectron case and have 
somewhat rounded edges. Examples of the muon energy 
spectra for some representative AKTW models display- 
ing these effects are shown in Fig. ^ for either beam po- 
larization configuration. Several things are to be noted 
in these Figures. The p, signals in the muon energy distri- 
bution vary over a wide range of both height and width 
depending upon the values of the p. and Xi masses and the 
production cross sections. The range of possible signal 
shapes relative to the background is varied, but gener- 
ally the signal is separable from the background in most 
models. We note again that the background is somewhat 
reduced compared to selectron production since there are 
fewer issues with beam remnants here. Again, we see that 
RH electron beam polarization leads to far smaller back- 
grounds than does LH beam polarization, as expected 
due to the diminished contribution from W-pair produc- 
tion. 

Of the 22 kinematically accessible models, 19(17) lead 
to signals with significance > 5 at these integrated lumi- 
nosities with RH(LH) electron beam polarization. Com- 
bining the LH and RH polarization channels, we find that 
19/22 models with accessible smuons lead to signals that 
meet our visibility criteria. We display a representative 
set of these models in Fig. ^for both beam polarizations. 
The three models that do not pass our discovery crite- 
ria have smuons with masses in excess of 241 GeV; this 
leads to a strong kinematic suppression in their cross sec- 
tions, and hence, very small signal rates. The S/B ratio 
is somewhat small for some of the models in the LH po- 
larization channel, as can be seen in Fig. |2^, and are not 
so easily visible at this integrated luminosity. However, 
they nonetheless pass our significance tests for discovery. 

It is again interesting to compare the AKTW models 
that have visible smuons at the ILC with the well studied 
case of SPSla'. Fig. ^ shows the muon energy spectrum 
we obtain after imposing our kinematic cuts in the case of 
SPSla' for both beam polarizations. As in our selectron 
analysis, we observe that the event yield for SPSla' is 
far larger than all the AKTW models we study here, in 
some cases by as much as a factor of order 50. Also, as in 
the previous analysis, two distinct shelves are observed 
since both jlL,R are being simultaneously produced. The 
muon energy distribution is slightly different from that 
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obtained for electrons in this model, not only because 
of the small differences in our cuts, but also due to the 
fact that the mixed final state /xl/z/j is not produced due 
to the absence of the i-channel contribution. Clearly, in 
comparison to the bulk of our models, it is rather trivial 
to discover and make precise determinations of the smuon 
properties in SPSla'. 

Note that 4 of the models with kinematically acces- 
sible smuons also have kinematically accessible lightest 
chargino states. However, since all of these charginos are 
rather close in mass to the LSP, i.e., within 5 GeV, the 
existence of the charginos does not constitute a large ad- 
ditional source of background and does not significantly 
affect the qualitative structure of the muon energy spec- 
tra. They could, however, modify the extracted values 
of the particle masses obtained from an analysis of the 
endpoints of the muon energy spectra and this possibility 
should be studied further. 

Interestingly, as in the selectron case above, a num- 
ber of models which do not have kinematically accessible 
smuons give rise to visible signals in the /i-pair analy- 
sis. This is just another example of the well-known phe- 
nomenon where SUSY is a background to itself. We find 
that there are 20(15) models which yield fake signals in 
the case of RH(LH) polarization. As in the previous anal- 
ysis, decays of other SUSY particles into muons, e.g., 
Xi W*Xi with W* "> ^v, can lead to the same ob- 
servable final state in both polarization channels. We 
present examples of such misleading signals in Fig. |2^. 
Note that in the case of RH beam polarization, the fake 
signature looks quite different than in a typical model 
which really has smuons present; there are no shelf-like 
structures and the muon energies are relatively low. This 
is to be expected when the final state muons are the 
result of a 3-(or more) body decay channel, and when 
chargino-neutralino mass splittings are small. However, 
for LH polarization, two representative fake models (la- 
belled 8324 and 39331 in the Figures) appear to mimic 
the smuon shelf-like feature. This is due to the fact that 
in these particular models, as will be discussed further 
below in Section 5.1, the W boson in the decay pro- 
cess is on-shell so that the final state muons are the result 
of true 2-body decays. 

In order to assist in the differentiation of models with 
real smuons from ones that do not, it is necessary to 
examine other kinematic distributions. Figures ^ and 
show the plji" distributions for the real smuon and fake 
models, respectively. Here we see that the models with 
real smuons generally lead to harder muons in the final 
state than do the counterfeit cases; this holds to some 
extent in the fake models where the charginos decay to 
on-shell W bosons. 
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FIG. 20: The muon energy distribution of the SM background after the smuon selection cuts described in the text have been 
imposed. The red(blue) solid curves correspond to the full background sample with 80% LH(RH) electron beam polarization, 
and the green(pink) dotted lines represent the contribution from the process e^e~ — > l~viviil' , I = h,t, with LH(RH) beam 
polarization, respectively. 250 fb~^ of integrated luminosity for each beam polarization has been assumed. 
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FIG. 21: Muon energy distribution: the number of events/2 GeV bin (combined signal and background) after imposing the cuts 
described in the text for several representative models, with RH(LH) beam polarization in the top(bottom) panel, assuming 
an integrated luminosity of 250 fb~^ for either polarization. The SM background is represented by the black histogram. 
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FIG. 22: Same as the previous Figure but now showing AKTW models that result in a small S/B ratio in the LH polarization 
channel in the //-pair analysis. 
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FIG. 23: Muon energy distribution: the number of events/2 GeV bin (combined signal and background) after imposing the cuts 
described in the text for the benchmark model SPSla', with RH(LH) beam polarization in the top(bottom) panel, assuming 
an integrated luminosity of 250 fb~^ for either polarization. The SM background is represented by the black histogram. 
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FIG. 24: Muon energy distribution: the number of events/2 GeV bin (combined signal and background) after imposing the 
cuts described in the text for representative models which lead to fake smuon signatures from chargino and neutralino decays. 
Here, we show RH(LH) beam polarization in the top(bottom) panel, assuming an integrated luminosity of 250 fb~^ for either 
polarization. The SM background is represented by the black histogram. 
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FIG. 25: The distribution for the same models shown in Fig. Ellwith real smuons. Here, we show RH(LH) beam polarization 
in the top(bottom) panel, assuming an integrated luminosity oi 250 fb~^ for either polarization. The SM background is 
represented by the black histogram. 
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FIG. 26: The p^^ distribution for the same models shown in Fig. t 
in the top(bottom) panel, assuming an integrated luminosity 
represented by the black histogram. 



I with fake smuons. Here, we show RH(LH) beam polarization 
[ 250 fb~^ for either polarization. The SM background is 
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Staus 

This analysis is similar to the other charged slepton 
analyses discussed above. We analyze the channel 

f+f- ^ T+r-x?X? , (10) 

that is, the signature is a tau pair plus missing energy. 
Staus are pair produced via s-channel 7 and Z exchange 
and receive no t-channel contribution. The left- and 
right-handed staus mix to form two mass eigenstates, 
which have mixing-dependent couplings to the Z boson. 
In contrast to the other charged slepton analyses above, 
the identification of the final state tau leptons is nontriv- 
ial, because the tau decays in the detector, predominantly 
into hadrons. 

We focus on the hadronic decays of taus into pions, 
T — > TTiyr', T — > pvr — !■ TT^TT^Vr', T — !■ iTTv-r, the latter be- 
ing a 3-prong jet, but also include the leptonic decays of 
the T. In the hadronic decay channel, taus are identified 
as jets with a charged multiplicity of 1 or 3, and with 
invariant mass less than some maximum value. Our tau 
selection criteria are as follows ^ (note that we em- 
ploy the notation tau-jet to describe the visible r decay 
products): 

1. We require 2 jets in the event, each with charged 
multiplicity of 1 (where the tau decays into a lep- 
ton, p, vr, or 37r-decay with 27r°s) or 3 (where the 
tau decays into 3 charged pions). 

2. The invariant mass of tau-jet, i.e., the visible tau 
decay products, must be < 1.8 GeV. 

3. If the tau-jet is 3-prong (charged multiplicity of 3), 
then none of the charged particles should be an 
electron or muon. 

4. If both tau-jets in the event are 1-prong, then we 
reject events where both jets are same flavor lep- 
tons, that is, an electron-positron or a muon pair. 
However we keep pairs of tau-jets that are, for ex- 
ample, an electron and a muon, or an electron and 
a pion, whereby a pion is defined as a charged track 
that is not identified as an electron or a muon. 

As an alternative analysis, we follow the above criteria 
and allow leptonic tau decays into muons, but reject taus 
that decay into electrons. This reduces contamination 
from photon-induced backgrounds. 

As mentioned above in our description of the SiD de- 
tector in Section , the current detector design does not 
allow for tracking, and hence does not have the capa- 
bility for particle ID, below 142 mrad. Thus muons at 
low angles are completely missed if they are too energetic 
to deposit energy into clusters. As we will see, certain 



7-induced processes constitute a significant background 
to stau production, particularly in the case where such 
energetic muons are produced but not reconstructed and 
the beam electron (or positron) receives a sufficient trans- 
verse kick to be detected. In this case, the final visible 
state is an electron and a muon, which would pass the 
standard tau ID preselection described above. The alter- 
native tau preselection criteria, which rejects the electron 
decay channel, eliminates this background at the price of 
reducing the signal correspondingly by roughly 30%. 

After these tau identification criteria are imposed, we 
employ cuts to reduce the SM background. Follow- 
ing |2^, we demand: 

1. No electromagnetic energy (or clusters) in the re- 
gion |cos6'| > 0.995. 

2. Two tau candidates as identified above, are 
weighted by their charge within the polar angle 
-0.75 < QrCosOr < 0.75. This reduces the Im- 
pair background. 

3. The acoplanarity angle must satisfy /S.(ff"^ > 40 de- 
grees. Here, since we demand two tau candidates, 
the acoplanarity angle is equivalent to tt minus the 
angle between the pr of the taus, A^"^"^ — n ~ 9t. 
The above requirement then translates to cos 9t > 
0.94. This cut reduces the W^-pair and 77-induced 
background. 

4. I cos 6'p„„^,„ J < 0.8. 

5. The transverse momentum of the ditau system be 
in the range O.OOS^s < Pt"' < 0.05^/i. This de- 
creases the 77-induced background. 

6. The transverse momentum of each of the tau can- 
didates be pt > 0.001-\/s. This cut is crucial to 
reduce the 77 and 67 background. 

7. The combined cut on J^P^^f ^-i^^ ^^^^i 
^p^^ < 0.001257^(1 + 5 sin A0^^) 

= 0.00125V^ (^1 + 5 y'l-cos^ 6*7^(11) 

is imposed. Here, X^P^f is the sum of the tau 
momenta projected onto the transverse thrust axis 
Ti, where the transverse thrust axis is given by 
the xy-components of the thrust axis. This last 
cut is necessary in the tau decay channel to further 
decrease the 77 background. 

As in the other slepton analyses, we histogram the re- 
sulting energy spectrum as well as PTvis in this case. 
We show the remaining SM backgrounds after these cuts 



38 



are imposed in Fig. |2^; the dominant background left af- 
ter the cuts stems ft'om 7-induced lepton-pair production 
processes. 

In Fig. ^ we also display the effect of the alternative 
tau ID criteria discussed above. This alternative tech- 
nique nearly completely eliminates the background since 
events where a beam electron/positron is falsely identi- 
fied as a tau decay product are rejected. Of course, as 
mentioned above this technique also reduces the signal 
by approximately 30%. Augmenting the detector with 
muon ID capabilities at lower angles could reduce the 
7-induced background without having to pay the price 
of introducing a restricted tau identification. As we will 
see below, a significant portion of the AKTW models 
have very low stau signal rates, and an improved track- 
ing capability could be crucial if in fact this portion of 
the SUSY parameter space is realized in nature. 

In 28 of the 242 AKTW models, the hghtest stau is 
kinematically accessible for pair production at the 500 
GeV ILC, and in one of these models the heavier stau 
partner can also be produced. The signal for stau pro- 
duction is somewhat different than in the case for selec- 
trons and smuons as the final state tau decays in the 
detector. In this case, we no longer have the distinctive 
shelf-like feature in the resulting energy spectrum of the 
reconstructed tau. The shape of these spectra is highly 
dependent on the mass difference between the stau and 
the LSP as shown in Fig. which displays the pure 
signal in 3 models before our selection cuts have been 
imposed. Here we see that small mass differences re- 
sult in a sharply peaked distribution at low tau energies, 
while a larger mass difference yields a flatter distribu- 
tion, albeit at a lower event rate. In principle, the search 
strategy, and hence the set of selection cuts, could be 
tailored to maximize the signal to background ratio once 
the stau— LSP mass difference is known. However, un- 
til the stau is discovered a general search strategy, such 
as that presented here, that applies for all mass regions 
must be employed. 

Of these 28 models, we find that 18 lead to signals 
which can be observed at the significance level S > 5. 
We also find that the heavier stau with ruf^ = 240 GeV 
is not produced with large enough event rates to be visi- 
ble above the SM background. In addition, in 3 of these 
28 models the mass difference between the lightest stau 
and the LSP is small enough such that the stau decays 
outside the detector, and it can be observed in our sta- 
ble charged particle search (described below in Section 
5.3). Of the 18 detectable stau models, 9(10) are visi- 
ble via our standard search criteria in the LH(RH) beam 
polarization configuration. The total (combining signal 
and background) energy distribution is shown for a 
representative set of these models for both beam polariza- 
tions in Fig. E9[ Here, we see that some of these models 



cleanly rise above the SM background, while others are 
just barely visible. The number of detectable stau sig- 
nals is greatly increased when we apply our alternative 
set of preselection criteria discussed above. One finds 
that 17(12) models are observable with LH(RH) electron 
beam polarization; the tau energy spectrum for a sam- 
pling of these models is displayed in Fig. In this case, 
we see that the signal is cleanly visible above the back- 
ground for all models. 

In all of the 18 models with observable staus, we find 
that the number of stau events that pass our cuts is dra- 
matically reduced compared to the event rate in Fig. |2^ 
before the cuts were applied. In addition, due to our 
cuts, the tau energy distribution is always peaked at low 
values, regardless of the stau-LSP mass difference, which 
ranges from 7—94 GeV in these 18 models. In the 7 
models where the signal is not observable, 2 are phase- 
space suppressed with mf^ > 240 GeV. The remaining 
5 models all have reasonable fi masses and the fi — Xi 
mass difference ranges from 42 — 108 GeV, but nonethe- 
less have a small production cross section due to stau 
mixing. 

Unlike many of the AKTW models we are examining, 
stau production at the ILC is straightforwardly observ- 
able in the case of the benchmark model SPSla'. This 
holds in either of the analysis channels as can be observed 
in Figs. |l] and Here we see that the stau signal is 
quite substantial and can be cleanly observed over the 
SM background for both choices of the electron beam 
polarization. 

We find that many AKTW models which do not con- 
tain kinematically accessible stau states nonetheless give 
rise to visible signatures with significance > 5 in this 
analysis, providing yet another example of SUSY being 
a background to itself. The tau energy distribution for a 
representative sampling of these SUSY background mod- 
els is presented in Fig. using the alternative set of 
kinematic cuts. We find that there are 29(28) models 
which yield fake signals with LH(RH) electron beam po- 
larization in our standard set of kinematic cuts. For our 
alternative analysis which rejects electrons in the final 
state, there are 30(28) models with false signatures for 
the LH(RH) polarization configuration. This analysis 
clearly has a very large number of false signals. We note 
that in every one of these "fake" models, xfxf^ X2X1J 
and X2X2 production is kinematically possible, and in one 
case selectron and smuon production is also viable. There 
are thus several sources of SUSY background which can 
lead to the same final state as that for stau pair produc- 
tion. 

In order to distinguish between stau production and 
these SUSY background sources, we investigate the vari- 
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FIG. 27: The energy spectrum of the SM background after the tau identification and stau selection cuts have been imposed. 
The solid blue and red lines are the full SM background for 80% right- (blue) and left-handed (red) electron polarization, 
the dashed lines (pink, right-handed, green left-handed) represent the dominant background source, e*7 ^ e^i/UjC^l^l' , I = 
e,jj,,T. Furthermore, the effect of eliminating the misidentification of beam electrons as tau decay products via the alternative 
tau ID described in the text is represented by the solid cyan line (80% right-handed beam polarization) and solid black line 
(80% left-handed polarization). 



able 

M^S = Er''+ E (12) 

1=1,2 

This variable is presented in Fig. ^ for both real and 
fake stau production. Here, we see that the false signals 
are slightly more peaked at lower values of M^q than 
does actual stau production. However, the distinction 
is not as clear as in the identification of selectron and 
smuon fake signatures discussed above, which makes use 
of the observable p^^* . This is because the full r energy 
is not carried by its visible decay products. We note 
that the p^*^ observable is not effective in distinguishing 
stau states from SUSY background sources, as in this 
case both the staus and the background sparticles have 
multi-body decays. 
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FIG. 28: Tau energy distribution stemming from stau pair production before our set of kinematic cuts have been imposed in 
three AKTW models. The red, blue, magenta histograms correspond to a stau-Xi mass difference of order 25, 100, and 0.5 
GeV, respectively. 
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FIG. 29: Tau energy distribution: the number of events/2 GeV bin after imposing the standard set of tau ID and f selection 
criteria described in the text for several representative AKTW models (signal and background combined). RH(LH) beam polar- 
ization is employed in the top(bottom) panel, assuming an integrated luminosity of 250 fb^^ for each polarization configuration. 
The SM background corresponds to the black histogram. 
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FIG. 30: Tau energy distribution: the number of events/2 GeV bin after imposing the alternative set of tan ID and f 
selection criteria described in the text for several representative AKTW models (signal and background combined). RH(LH) 
beam polarization is employed in the top(bottom) panel, assuming an integrated luminosity of 250 fb~^ for each polarization 
configuration. The SM background corresponds to the black histogram. 
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FIG. 31: Same as in Fig. ^ except now for SPSla'. 
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FIG. 32: Same as in Fig. |^ except now for SPSla'. 
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FIG. 33: Tau energy distribution: the number of events/2 GeV bin after imposing the alternative set of tau ID and f 
selection criteria described in the text for several representative AKTW models that give a fake signal in this channel. RH(LH) 
beam polarization is employed in the top(bottom) panel, assuming an integrated luminosity of 250 fb~^ for each polarization 
configuration. The SM background corresponds to the black histogram. 
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FIG. 34: EfTective mass (as employed at the LHC) of the products from the decays of objects identified as tau pairs produced 
from stau production (top panel) and from other, non-stau SUSY sparticles (bottom panel). Here we present the number of 
events/2 GeV bin assuming RH polarization and 500 fb~^ of integrated luminosity for several models shown in comparison to 
the SM background represented as the black histogram. 
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Sneutrino Pair Production 

We now examine the neutral slepton sector, i.e., sneu- 
trinos, which provides another potential handle for dis- 
tinguishing between models. For all three sneutrino fam- 
ilies there is the usual Z boson exchange contribution 
to the production amplitude in the s-channel, while for 
electron sneutrinos there is an additional ^-channel graph 
due to xf 2 exchange. If the charginos are heavy, then 
the Z-exchange graph dominates for all three generations 
and the resulting production cross section is determined 
solely by the amount of available phase space. 11/242 of 
our AKTW models have electron or muon sneutrino pairs 
which are kinematically accessible at -/s = 500 GeV, 
while 18/242 models contain accessible tau sneutrinos. 
In one of the models the sneutrino is the LSP. 

Sneutrinos, being neutral and weakly-interacting, are 
essentially only visible through their decay modes, of 
which there are several possible channels to consider: 
(i) vi viXi l^-rgely dominates in most cases, but leads 
to an invisible final state which, by itself, is clearly useless 
for either discovery or model comparison, (ii) i>i — > W£ 
is kinematically forbidden as an on-shell mode when 
i — e, fj, in all of our AKTW models and thus the corre- 
sponding 3-body branching fraction mediated by off-shell 
W bosons is very small. However, due to large f mixing, 
this 2-body mode is allowed for 6 of our models in the 
case of £ = r. When both the W and r decay hadron- 
ically, we can search for final states with multiple jets 
plus missing energy in this case, (iii) vg — > can also 
occur, with the subsequent decay X2 ~^ XiJj via a Z or 
Higgs boson. This occurs in 1(3) models in our sample 
when £ = ej ^.{t). However, in this case the resulting jets 
are likely to be relatively soft, due to a smaller — 
mass difference, making this mode difficult to observe 
above background, {iv) vi —* Xi^^ is accessible in 1(6) 
of these models when £ = e/nir), and leads to a final 
state of multiple jets plus two charged leptons plus miss- 
ing energy. As before, it is probable that these jets will be 
soft due to a smaller mass splitting between the chargino 
and the LSP and will most likely be difficult to observe 
depending upon the details of the rest of the spectrum. 

We first study the final state + missing en- 

ergy. This final state results from the decays 

u ^ iy[->jj7x?, 

i> ^ IXi^ Ijjxl ■ (13) 

For our signal selection, we require: 

1. Precisely one opposite sign lepton pair and two jet- 
pairs and no other charged particles in the event. 

2. No particles/clusters below the angular region of 
100 mrad. 



3. Missing energy to be > 2m^[)^jj^. We take 
m^o = 46 GeV, which is the current (weak, yet 
model-dependent) bound on the mass of the light- 
est neutralino ||^ . This bound arises from the invis- 
ible decay width of the Z boson and holds unless 
the Xi is very fine-tuned to be a pure Bino state 
and thus has no couplings to the Z Q. However, 
in order to estimate the effect on the background if 
this bound is increased, we perform a second anal- 
ysis with TOj^o = 100 GeV. 

4. In order to eliminate background that originates 
from very soft leptons or jets, we demand -Ejct,; > 
O.OlVs- 

These cuts effectively remove most of the SM back- 
ground as can be seen in Fig. |3^, which displays the miss- 
ing energy distribution for the background sample. Here, 
we see that at large values of missing energy, the domi- 
nant background remaining after the cuts arises from the 
process 77 cc, bb. Unfortunately, the sneutrino signal 
rates are also small. Fig. ^ presents the results of our 
analysis: none of the sneutrino models rise above the 
background, but several of the 'fake' models, where the 
signals arise from other SUSY particles, do appear. 4(3) 
fake models yield visible signals in the case of RH(LH) 
beam polarization. The counterfeit signals here are due 
to chargino and neutralino production and decay. We 
find that increasing the minimum LSP mass to 100 GeV 
docs not improve these results. 

We also study a second channel, with 6 jets and miss- 
ing energy in the final state. This is produced from the 
decay i) — * Wf — > jjjXi- The cuts and observables are 
similar to those of the 4j2£+ missing energy analysis, 
with the obvious substitution that we demand precisely 
6 jets and no other charged particles to appear in the 
event. We find that there is little to no SM background 
in this channel. However, we also find that none of our 
models are observable in this channel. 

An additional possible way to observe sneutrinos is via 
the radiative process e^e^ — > C/eiye + j. This may be 
particularly useful in the case where the decay channel 
j>f — > i/'f -|- Xi dominates. This reaction leads to a fi- 
nal state with a photon and missing energy and is thus 
similar to radiative LSP pair production, which we will 
discuss in detail below in Section 6.2. We find that radia- 
tive sneutrino production generally occurs with a smaller 
cross section than its LSP counterpart. As we will see be- 
low, the SM background from e^e~ —^ vD"f are generally 
too large to see this radiative process. 

In the case of the SPSla' benchmark model, the sneu- 
trino pair production mode is invisible as the sneutrinos 
dominantly decay into the uxi final state, as in most of 
our models here. 

Taking these results together for these various sneu- 
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FIG. 35: Missing energy distribution in 2 GeV bins for the SM background after the sneutrino selection cuts for the 4jet +21 
channel have been imposed. The blue(red) histogram corresponds to 80% RH(LH) electron beam polarization. The green 
dotted curve corresponds to the dominant background source, 77 cc,bb. 250 fb~^ of integrated luminosity was assumed for 
each polarization channel at 500 GeV. 



trino analyses, we conclude that the direct observation 
of i>i production is very difficult, if not essentially hope- 
less for the set of AKTW models. 
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FIG. 36: Missing energy distribution for the sneutrino 4j+lepton pair analysis: the number of events/2 GeV bin for several 
the fake models, with RH(LH) beam polarization in the top(bottom) panel, assuming an integrated luminosity of 250 fb~^ for 
either polarization. The SM background is shown as the black histogram. 
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CHARGING PRGDUGTIGN 

The chargino sector is simpler than that of the ncutrali- 
nos as in the CP-conserving MSSM; it depends on only 
three real Lagrangian parameters at tree- level: M2,/i, 
and tan/3. The resulting mass eigenstates, x?2) 



thus general admixtures of the charged Wino and Hig- 
gsino weak states. Figure 37 displays the Wino/Higgsino 
content of the lightest chargino, Xi^j in the 53 AKTW 
models that contain kinematically accessible charginos at 
y/s = 500 GeV. We see that the lightest chargino tends 
to be an almost pure Wino or Higgsino state in most of 
our models. In e+e~ collisions, such particles can be pro- 
duced via two mechanisms; s-channel 7, Z exchange pro- 
duces either pure Wino or Higgsino pairs but no mixed 
Wino-Higgsino final states as the analogous WZH^ cou- 
pling is absent. In addition, the t— channel sneutrino ex- 
change amplitude produces pairs of charged winos only. 
Clearly, the cross section for chargino pair production 
not only depends on the eigenstate masses but also on 
the various mixing angles present in the chargino sector. 
At ^/s = 500 GeV, pairs are typically too heavy to 
pair produce, so we will consider only xt P^i^' production 
in our analysis below. p9[| 

Once produced, the detailed nature of the xt de- 
cays critically depends upon the mass difference Am^ — 
m-± — m^o, with the latter being the LSP. In all cases, 

xt decay can proceed via either a W , Xi ~^ VVxi, or 
through an intermediate slepton, e.g., xt ~^ £v,i' — > 
i^Xi- seen above in Fig. ^, and in Figs. |3^ and ^ the 
mass spectrum of the models that have chargino states 
kinematically accessible at ^/s = 500 GeV is such that 
the distribution of values for Am^ are concentrated in 
the region < 5 GeV. The variation in Am^ yields several 
distinct signatures for xt production and thus all possi- 
ble values of Am-^ must be considered when performing 
our analysis. For example, if Am^ > Myy then on-shell 
W bosons can be produced and may be identified either 
through their leptonic or hadronic decay modes. We thus 
consider channels such as xt ~^ j.'J ^misu, with the dijets 
reconstructing to the W mass, or xt ~^ M-E-miBs, with 
the latter mode also covering decays through the slep- 
ton channels: xt ~* fl^v^fi^iy —> fiE^^^^. Hence, in this 



region, we search for the final states Xi^Xi 



2j fs^ -Emiss or fJ,~^iJ,~ ii'miss- For smaller values of Am, but 
still greater than a few GeV, we search for the same fi- 
nal states although the dijets will no longer reconstruct 
to Mw- A more difficult region is reached when Am^ 
is only a few GeV or less, as then the visible part of 
the xt decays are very r-like. Fig. ^ shows this region 
of chargino-LSP mass splitting for the range Am-^^ < 6 
GeV as a function of the chargino mass; we see that this 
region comprises the bulk of our models with accessible 
xt states at -^i = 500 GeV. 



At the other end of the spectrum, we must consider 
the case of small values of Am^ <, 1 GeV. The branching 
fractions for the xt decay channels in this case are pre- 
sented in Fig. If Am^ is very small, 100 MeV, then 
the chargino is long lived, and will travel many meters 
before decaying into an electron and missing energy. In 
such a case, we perform a massive stable charged particle 
search, determining the velocity of the xt momen- 
tum and energy measurements. As Am^ increases from 
this tiny value and the thresholds for xt decay into the 
fi and pion(s) are passed, the chargino lifetime decreases 
substantially and the chargino now decays to soft charged 
particles. In this mass range there are two possible search 
techniques: one can either look for decays in the detector 
from (semi-) long-lived xt states, or tag these soft decays 
via photon emission off of the initial and final state par- 
ticles. The latter corresponds to the radiative process 
e+e^ XiXil nil 07 f-iid is the approach we will 
pursue below. 

At = 500 GeV, we find that 53/242 of the AKTW 
models have kinematically accessible charginos. Fig- 
ure ^ shows that in all but two cases (which we label 
here as models 8324 and 39331) the models populate the 
region Am^ < 5.5 GeV; for the two exceptions we see 
that Arrij^ > 100 GeV. Interestingly, we note that mod- 
els with small Am < 1 GeV tend to have a large Wino 
content, while those in the range Am^ ~ 4 '--^ 5 GeV are 
found to have a large Higgsino content as can be seen in 
Fig. |§. 

We now discuss our analyses for each region of Am^^. 



Non-Close Mass Case 



We first examine the case where 



Am-j;; 



m^i 



m^o > 1 GeV. 



(14) 



As already mentioned above, there are several possible 
final states that can be studied in this mass region. Anal- 
ysis techniques have been developed for the cases where 
the chargino decay proceeds through on- or off-shell W 
bosons (with W — > jj or W ^ Ivi) or via sleptons {e.g., 
p, — > /iXi)- We discuss each of these in turn. 



Chargino Decays via On-Shell W bosons 

This analysis applies to the case Am^ > Mw- Here, 
we examine four-jet final states, stemming from the de- 
cays of the chargino pair into W bosons with subsequent 
decays into quark pairs. 



w 



-w-xlxl 



(15) 
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Chargino Mixing 




iuiir 

FIG. 37: Wino/Higgsino content of the xt ^'^^ 53 models that have kinematically accessible charginos at ^/s = 500 GeV. 
U{V) is the left(right) diagonalizing matrix. In the lower left (upper right) corner of the Figure, the physical chargino is 
dominantly Higgsino (Wino). 
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FIG. 38: Lightest chargino-LSP mass difference for the region < Am^ < 800 GeV. Note that the chargino states in the 
models where Amx > 150 GeV are not kinematically accessible at sqrts — 500 GeV. 
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FIG. 39: Lightest chargino-LSP mass diflerence for the close mass region < Am^ < 1 GeV. 
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with 



qq. 



(16) 



Taking the hadronic decay mode of both W bosons yields 
the final state with the largest statistical sample. 

As always, the SM background is significant. In order 
to reduce the background, we demand, as expanded and 
adapted from @ ||: 

1. There be precisely 4 jets in the final state and no 
other charged particles. As mentioned in section , 
we employ the JADE jet algorithm in the E scheme, 
with ?/cut = 0.05. This choice of ycut avoids the 
situation where soft gluons produce too many jets. 

2. E^i^ > in the backward direction. Here, the back- 
ward direction is defined with respect to the thrust 
axis, and corresponds to the hemisphere with the 
lesser amount of energy. This cut is designed to 
reduce SM background from Z pair production, 
where one Z boson decays into neutrinos, and the 
other Z decays into quarks, which then radiate hard 
gluons. 

3. In the forward direction, the visible energy is con- 
strained to be E^i^ < ■^y/s — n^x'i,'am- the 
case of stau production, we take m^o^^-^^ — 46 GeV, 
which is the current bound on the mass of the light- 
est neutralino Q in the case where the Xi is not a 
pure Bino eigenstate. 

4. The visible energy is constrained to be E^^^ < 1 
GeV in the region 0.9 < | cos6'| < 0.99. This is to 
decrease the W pair background which is strongly 
peaked in the forward direction. 

5. The acoplanarity angle satisfy Af?!)^'''^'''' j'^'p"'' > 30 
degrees. Since we demand that two jet pairs re- 
combine into W bosons, the acoplanarity angle is 
equivalent to tt minus the angle between the pt of 
the W bosons, i.e., Acj)'"''"''' = tt-Ot. This sig- 
nificantly reduces the W pair and 77 background, 
since the W bosons from the chargino decays are 
accompanied by missing energy from the LSP. 

The standard search analysis for this final state is 
based on the energy distribution of the jet-pairs which 
reconstruct into a W boson. As in the case of the se- 
lectron analysis, this distribution should have a box- like 
shape with a shelf and sharp endpoints in the presence 
of a high statistics, background-free, perfect detector en- 
vironment with the absence of radiative effects. Here, 
the 2-body decay is taken to be ^ + X? and the 
expressions in Eqs. (0)-(||) need to be adapted to include 



the massive W boson. In this case, one can solve the 
equations for the chargino and LSP masses and finds 



2 A±B 

Tfi A- 

xt 2{E^,^ + E^in)ys ' 



(17) 



rn^± + ^ 1 - 2 -= ) , (18) 



where 



Ern^-vE^ 



maxJ^min i 



B = 



£■2. ) 



(19) 
(20) 



and where -Emax, ^'min arc determined experimentally. 

We refrain from presenting the remaining SM back- 
ground after these cuts are imposed, as only a handful of 
events pass the cuts. As mentioned above, only 2 models 
in our sample lie in the kinematic region Am-^ > Mw- 
The jet-pair energy distribution for these two cases is 
displayed in Fig. |4^ for left- and right-handed electron 
beam polarization. Here, we see that the overall event 
rate that survives the kinematic cuts is not large, but the 
signal cleanly towers above the even smaller background. 
We can see the effects of the cuts and the detector en- 
vironment in these cases, as the shape of the spectrum 
does not display the shelf-like behavior discussed above. 
Note that an additional model (labelled 1822), which has 



At 



1 GeV, also passes the kinematic cuts for this 



analysis. This model yields a smaller event rate than the 
cases with on-shell W bosons, but populates a different 
region of the spectrum. However, this model contains 
a light : Xi 2 3 sector and the signal that passes our 
cuts here is due to the production of these heavier gaug- 
ino states and not from the Xi and is thus a fake. Three 
out of 53 AKTW models with accessible charginos are 
thus visible with a significance 5 > 5 in this channel. 
We note that fake signals from the production of other 
SUSY particles do not satisfy our visibility criteria in 
this channel, except for model 1822. It would thus seem 
that this analysis is relatively free from Supersymmetric 
backgrounds, at least in the case of the AKTW models. 

In summary, if Am^ is large enough to produce on- 
shell W bosons, this is clearly a very clean channel. 



Chargino Decays via Off-Shell Ws and/or Sleptons 

We search for three final states in the kinematic region 
Mw > Amy > 1 GeV where the W boson is produced 
off-shell in the decay: four jet events plus missing 
energy, two jets and a lepton plus missing energy or two 
leptons plus missing energy. In order to avoid the large 
SM background from the beam remnants in 76^, 77 re- 
actions, we require that the final state leptons be muons. 
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FIG. 41: Jet-pair energy distribution: the number of events/2 GeV bin after imposing the full set of cuts discussed in the text 
for chargino production and on-shell decays to a boson for the three models which are visible in this channel. RH(LH) beam 
polarization is employed in the top(bottom) panel, assuming an integrated luminosity of 250 fb~^ for either polarization. The 
SM background is shown as the black histogram. 
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In the fully leptonic decay channel, the kinematic cuts 
we employ to distinguish signal from background are 
very similar to the slepton searches described in Sec- 
tion above. As in the case of smuon production, one 
searches for a structure above the SM background in the 
final state muon energy distribution. Here, however, the 
signal distribution is not expected to display the by-now 
familiar shelf-like behavior, due to the 3-body nature of 
the chargino decay. Examples of the muon energy spec- 
tra for some representative AKTW models are shown in 
Fig. ^ for both beam polarization states. Here, we see 
that the Efj_ spectrum for the signal varies greatly in size 
and shape between the various models depending on the 
value of Arrijj- and the production cross section, but is 
nonetheless clearly separable from the SM background 
for these cases. Comparing with Fig. |T]we note that the 
SM background has a similar shape, but is slightly larger 
throughout the spectrum in this analysis. As usual, RH 
electron beam polarization leads to a smaller SM back- 
ground since the t-channel contribution to W pair pro- 
duction is suppressed in this case. 

Of the 53 models with kinematically accessible xf 
states, we find that 12(11) lead to visible signals over 
the background at a significance of 5 > 5 for RH(LH) 
electron beam polarization. Combining the two polariza- 
tion channels, a total of 14/53 models meet our visibility 
criteria. From Fig. ^ we note, however, that S/B can 
be small enough in some cases to render a detailed study 
of the chargino properties difficult. 

SUSY can be a substantial background to itself in 
this channel, with smuon production being a particularly 
large background source. We find that 14(12) models 
yield fake signals that pass our visibility criteria for a 
RH(LH) polarization configuration. We note that in all 
cases, the counterfeit signal indeed arises from the pro- 
duction of smuons. The muon energy distribution for 
some representative examples of such misleading signals 
are presented in Fig. From the Figure, we see that 
the signal tends to have a shelf-like behavior, as would 
be expected for smuon production, and thus looks quite 
different than the case of chargino production. 

We now turn to the fully hadronic channel, where the 
final state is 4 jets plus missing energy. We employ the 
following kinematic cuts (based on, e.g., ||35||): 

1. There be precisely 4 jets in the final state and no 
other charged particles. 

2. No tracks (or clusters) be present below an angle 
of 100 mrad. This reduces photon-initiated back- 
grounds. 

3. Missing energy is constrained to be > 0.5^/s. This 
favors the signal, which contains a large amount 
of missing energy compared to many background 
sources. 



4. In the forward direction, the visible energy is con- 
strained to be -EviB < — TOxJ.min- ^g^i'^i 
take m^o = 46 GeV. However, in order to es- 
timate the effect on the background if this bound 
is increased, for example by future studies at the 
LHC, we perform a second analysis with m^o^^^^^ — 
100 GeV. 

5. We require precisely two jets in each hemisphere as 
determined by the thrust axis. This cut eliminates 
jets stemming from t decays arising from tau pair 
production, where one t has a one-prong decay, and 
the other is 3-prong. 

6. We reconstruct the off-shell W bosons by coalescing 
the 4 jets into 2, one for each W boson. We force 
this by adjusting the j/cut parameter of the JADE 
jet finding algorithm until two jets are found (see 
Section for details of the algorithm), and then re- 
quire the resulting dijet invariant masses to be > 2 
GeV. 

The first observable we consider for this channel is the 
missing energy distribution, where we expect a peak from 
the signal at large values of i?miss- This spectrum is pre- 
sented in Fig. |44| where the black histogram corresponds 
to the SM background, as usual. We find that a partic- 
ularly troublesome background arises from the process 
77 ~^ QQ- It is clear from the Figure that this back- 
ground reaction stubbornly yields a significant event rate 
even after the above cuts are imposed. The missing en- 
ergy spectra for representative AKTW models are also 
shown in the Figure for both polarization states and by 
eye seem barely visible above the background. However, 
the statistical sample is large, and nonetheless, chargino 
production in several AKTW models are observable in 
this channel. We find that 9(31) models are observable 
with significance > 5 in the RH(LH) polarization state. 
All models that are visible with RH beam polarization 
are also observed with LH beam polarization. If we in- 
crease the minimal value of the bound on the LSP mass 
to 100 GeV as mentioned in the kinematic cuts above, 
we find that 4 additional models satisfy our visibility cri- 
teria. We note that the production of other Supersym- 
metric particles do not pass our cuts in this channel and 
hence there are no 'fake' signals. 

We next examine the distribution of the two dijet in- 
variant masses to see if the photon-induced background 
is less problematic for this observable. Figure ^ displays 
the SM background for this case for both electron beam 
polarizations. Here, we see that the background from 
two-photon initiated processes is still significant, with the 
dominant channel passing the cuts being 77 — > qq. The 
results for the AKTW models are presented in Fig. ^ for 
both electron beam polarizations. We find that with RH 



56 



E(mu+,mu-)— chargino 



RS324.aida 
R10197.aida 
R13960.aida 

- R14343.aida 

- R20992.aida 

- R27712.aida 

- mixRI.aida 



RS324.aida 
Entries : 3954 
Mean : 94.722 
Rms : 64.785 

R10197.aida 
Entries : 7090 
Mean : 60.480 
Rms : 60.983 



R1 3960.aida 
Entries : 4608 
Mean : 82.401 
Rms : 65.727 




R1 4343.aida 
Entries : 3312 
Mean : 101.06 
Rms : 68.353 



R20992.aida 
Entries : 3482 
Mean : 98.083 
Rms : 68.196 



R27712.aida 
Entries : 3230 
Mean : 104.07 
Rms : 67.298 

mixRI .aida 
Entries : 31 46 
Mean : 105.40 
Rms : — 67.367 



E(mu+,mu-)— chargino 



- L8324.aida 

- L8742.aida 

- L10197.aida 
LI 2843. aida 
L13960.aida 

- L39331 .aida 

- mixLl.aida 



L8324.aida 
Entries : 18665 
Mean : 105.84 
Rms : 66.153 



L8742.aida 
Entries : 16919 
Mean : 1 1 1 .26 
Rms : 66.969 




L10197.aida 
Entries : 18501 
Mean : 103.70 
Rms : 68.763 

L1 2843.aida 
Entries : 16713 
Mean : 112.23 
Rms : 66.772 



FIG. 42: Muon energy distribution: the number of events/2 GeV bin after imposing the full set of cuts discussed in the text 
for chargino production for representative models which are visible in this channel. RH(LH) beam polarization is employed in 
the top(bottom) panel, assuming an integrated luminosity of 250 fb~^ for either polarization. The SM background is shown as 
the black histogram. 
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FIG. 43: Muon energy distribution: the number of events/2 GeV bin after imposing the full set of cuts discussed in 
text for representative models which fake a chargino signal in this channel. RH(LH) beam polarization is employed in 
top(bottom) panel, assuming an integrated luminosity of 250 fb~^ for either polarization. The SM background is shown as 
black histogram. 
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FIG. 44: Missing energy distribution: the number of events/2 GeV bin after imposing the fuh set of cuts discussed in the 
text for the fully hadronic chargino decay channel for representative models which are visible in this channel. RH(LH) beam 
polarization is employed in the top(bottom) panel, assuming an integrated luminosity of 250 fb^^ for either polarization. The 
SM background is shown as the black histogram. 
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polarization none of the AKTW models with chargino 
decays into off-shell W bosons are visible over the back- 
ground in this observable. In the case of LH electron 
beam polarization, we find that 2 such models (labelled 
as 12843 and 14343) are detectable with a significance 
> 5. These models yield an excess of events in the first 
two bins of the distribution; this excess is not visible by 
eye, but is statistically significant due to the large sample 
size. The two models with chargino decays into on-shell 
W bosons (labelled as 8324 and 39331) display a clear 
signal for both beam polarizations as shown in the Fig- 
ure. We see that the invariant Mjj spectrum is broader 
for these two models and yields a high event rate at large 
invariant masses. If the minimum value of the LSP mass 
is raised to 100 GeV in our kinematic cuts, we find that 
only the on-shell W boson decays are visible above the 
background. One additional model (labelled as 1822) is 
also distinguishable from the background. In this case, 
however, it is due to the pair production of X2 states, 
with their subsequent decays into Z + ^ 2jets -|- x?, 
that passes our cuts and provides a fake signal. 

Further attempts to decrease the SM background in 
the 4-jet channel prove to be difhcult and tend to remove 
the signal as well as the background. This is because the 
characteristics of the signal and remaining background 
are similar for variables such as the missing energy spec- 
trum, acoplanarity, and pTvia (see Fig. In partic- 
ular, both signal and background distributions peak at 
low values of pTvis and acoplanarity and at high values of 
missing energy as discussed above. Previous searches for 
charginos in the literature {e.g., |35 ) have employed ad- 
ditional cuts on PTvis and/or acoplanarity. In particular. 



we find that an additional restriction on the transverse 
momentum, 

7. PT„. > o.oeyi 

effectively reduces the background as shown in Fig. 
but also removes the signal for all of the AKTW mod- 
els with chargino decays into off-shell W bosons. We 
show the effect of this additional cut on transverse mo- 
mentum in the missing energy and jetpair invariant mass 
distributions in Figs. ^ and ^ for both electron beam 
polarizations. Here, we display all of the AKTW models 
which yield a visible signal with significance > 5. We see 
that the signal for the models where the chargino decays 
into on-shell W bosons (labelled as 8324 and 39331) tow- 
ers above the background for both observables with both 
beam polarizations. As we saw above, the visible sig- 
nal for model 1822 is due to the production of X2 states 
and is thus fake. We emphasize that none of the AKTW 
models with chargino decays into off-shell W bosons are 
observable once this additional pt cut is applied. 
Lastly, we examine the mixed decay channel. 



which can proceed via (if kinematics allow) 

Xtxl ^W*^Xl+^^^'>^.,W*^Xl+^^^^^. (22) 

with 

W*^^qq, (23) 



xtxi 



W*'^W* 



(24) 



where one of the virtual W bosons decays hadronically 
while the other decays leptonically into a muon. 
For this channel, we employ the cuts 

1. There be 2 jets (with no muonic component) plus 
one muon with no other visible particles in the final 
state. 

2. There be no tracks or clusters of energy within 100 
mrad of the beampipe as the signal is peaked at 
wide angles. 

3. The visible energy satisfy E^^^ < ^^/s — rrij^o^^^^ in 
the forward direction. As above, we take TOco = 

. However, in order to estimate the effect 
on the background if this bound is increased, we 
perform a second analysis with m^o = 100 GeV. 



Xi Xi 



mxi 



f^'^iJ.Xi , 



4. The invariant mass of the jet-pair be larger than 2.4 
GeV. This cut is used to eliminate jets stemming 
from T decays. 

Here, we examine the energy and invariant mass of the 
jet-pair. We find that the same issues discussed above 
in the 4-jct channel regarding potential additional cuts 
based on acoplanarity and transverse momentum are also 
relevant in this case and thus these cuts are not em- 
ployed in our analysis. The background remaining after 
our cuts are imposed is presented in Fig. ^ for the in- 
variant mass distribution, where the dominant remaining 
SM background processes are 77 — > qq, 77 t^t" and 
e^e~ — > qq'Wi- We see that the background distribution 
is roughly the same for both beam polarizations as is to 
be expected for 77 induced processes. 

The invariant mass spectrum for the case where the 
charginos decay to off-shell W bosons is displayed in 
Fig. ^ for both beam polarizations. In both cases, the 
signal rises above the background in the region of smaller 
(< 60 GeV) invariant masses of the jet pair, as is ex- 
pected due to the off-shell nature of the W bosons. The 
model labelled as 1822 also shows a visible signature, 
with a peak located at Mjj ~ 80 — 90 GeV. As we saw 
above, this is due to X2 production in this model with 
the subsequent decay X2 ~^ ^ ~^ Xi with the Z decaying 
(21) hadronically and is a false signal. The Mjj distribution 
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for the 2 AKTW models where the charginos decay to 
on-shell W bosons is shown in Fig. Here, we would 
expect to see a peak above the SM background in the 
distribution around Mw, and indeed, that is the case. 
In summary, we find that 23(35) of the AKTW models 
with kinematically accessible charginos lead to signals in 
this observable with a visibility significance 5 > 5 for 
RH(LH) electron beam polarization at these integrated 
luminosities. We note that none of the AKTW models 
are visible over the background in the case where the 
minimum value of the LSP mass is increased to 100 GeV 
as described in our cuts. 

The second observable we examine in this analysis is 
the energy of the jet pair which is displayed in Fig. |5j 
for several AKTW models where the chargino decays 
to an off-shell W boson. Again, we see that the sig- 
nal rises above the background for lower values of Ejj 
(<^ Mw), except for the case of model 1822 which is a 
fake as described above. For comparison, the results for 
the two models which have decays to on-shell W bosons 
are shown in Fig. |5^, where we see that the Mjj spectrum 
is peaked at larger values in this case. For this observ- 
able, we find that 26(35) of the AKTW models meet our 
visibility criterion. We note that two more models are 
visible with RH polarization in this observable compared 
to the Mjj distribution discussed above. Again, none of 
the models are visible when the minimum value of the 
LSP mass is increased in the analysis. 

Except for model 1822, we note that there are no fake 
signals from the production of other SUSY particles for 
either observable in this channel. 

We now compare these results to those for the case 
of the well-studied benchmark point SPSla'. Figures |5^ 
and 57 display the jet pair invariant mass spectrum and 
energy distribution, respectively, for both polarization 
choices. The chargino in this model decays to an on- 
shell W boson and has a large production cross section; 
both of these features are observable in the Figures. The 
signal for this model is clearly visible above the SM back- 
ground and there is a peak at Mjj ~ Mw in the invariant 
mass distribution. 



61 




FIG. 45: Distribution of dijet invariant masses from the remaining SM background after the chargino 4-jet selection cuts listed 
in the text have been imposed. This is generated for 250 fb~^ of SM events with 80% right-handed (solid blue line) and 80% 
left-handed (solid red line) electron beam polarization, and unpolarized positron beam at ^/s = 500 GeV. The dashed green 
line shows the main processes contributing to the background, 77 — > qq, which is independent of the beam polarization. 
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FIG. 46: Jet-pair invariant mass distribution: the number of events/2 GeV bin after imposing the full set of cuts discussed in 
the text for the fully hadronic decays of the chargino for representative models which are visible in this channel. RH(LH) beam 
polarization is employed in the top(bottom) panel, assuming an integrated luminosity of 250 fb~^ for either polarization. The 
SM background is shown as the black histogram. 
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FIG. 47: Transverse momentum distribution of remaining SM background after the chargino 4-jet selection cuts listed in the 
text have been imposed. This is generated from 250 fb~^ of SM events with 80% left-handed electron beam polarization, and 
unpolarized positron beam at = 500 GeV (solid blue line). The dashed purple and red lines show signal events produced in 
two AKTW models that are representative for the class of models with Arrijj of the order of a few GeV. 
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FIG. 48: Distribution of dijet invariant masses from the remaining SM background after the chargino 4-jet selection cuts listed 
in the text have been imposed with the additional cut on transverse momentum. This is generated from 250 fb~^ of SM events 
with 80% right-handed (dashed pink line) and 80% left-handed (dashed green line) electron beam polarization, and unpolarized 
positron beam at ^/s = 500 GeV. The solid lines are as in fig. Ha 
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FIG. 49: Missing energy distribution: the number of events/2 GeV bin after imposing the fuU set of cuts discussed in the text 
for the fully hadronic chargino decay channel, including an additional cut on transverse momentum, for representative models 
which are visible in this channel. RH(LH) beam polarization is employed in the top(bottom) panel, assuming an integrated 
luminosity of 250 fb~^ for either polarization. The SM background is shown as the black histogram. 
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FIG. 50: Jet-pair invariant mass distribution: the number of events/2 GeV bin after imposing the full set of cuts discussed in 
the text for the fully hadronic chargino decay channel, including an additional cut on transverse momentum, for representative 
models which are visible in this channel. RH(LH) beam polarization is employed in the top(bottom) panel, assuming an 
integrated luminosity of 250 fb~^ for either polarization. The SM background is shown as the black histogram. 
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FIG. 51: Distribution of dijet invariant mass from the remaining SM background after the chargino 2-jet plus muon selection 
cuts listed in the text have been imposed. This is generated from 250 fb~^ of SM events with 80% right-handed (solid blue line) 
and 80% left-handed (solid red line) electron beam polarization, and unpolarized positron beam at = 500 GeV. The other 
dashed, and solid lines show the main processes contributing to the background, 77 qq (dashed green line), 77 t^t~ 
(dashed pink line), which are independent of beam polarization, and e^e'^ — > qq'Wi, for 80% right-handed (solid cyan line) and 
80% left-handed (solid black line) electron polarization. 
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FIG. 52: Jet-pair invariant mass distribution: the number of events/2 GeV bin after imposing the full set of cuts discussed in 
the text for the 2-jet +// channel with Am<i- < Mw for representative models which are visible in this channel. RH(LH) beam 



polarization is employed in the top(bottom) panel, assuming an integrated luminosity of 250 fb 
SM background is shown as the black histogram. 
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FIG. 53: Jet-pair invariant mass distribution: the number of events/2 GeV bin after imposing the full set of cuts discussed in 
the text for the 2-jet +// channel with Am<i- > Mw for representative models which are visible in this channel. RH(LH) beam 
polarization is employed in the top(bottom) panel, assuming an integrated luminosity of 250 fb~^ for either polarization. The 
SM background is shown as the black histogram. 
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FIG. 54: Jet-pair energy distribution: the number of events/2 GeV bin after imposing the full set of cuts discussed in the 
text for the 2-jet +fi channel with Am^ < Mw for representative models which are visible in this channel. RH(LH) beam 
polarization is employed in the top(bottom) panel, assuming an integrated luminosity of 250 fb~^ for either polarization. The 
SM background is shown as the black histogram. 
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FIG. 55: Jet-pair energy distribution: the number of events/2 GeV bin after imposing the full set of cuts discussed in the 
text for the 2-jet channel with Am^ < Mw for representative models which are visible in this channel. RH(LH) beam 



polarization is employed in the top(bottom) panel, assuming an integrated luminosity of 250 fb 
SM background is shown as the black histogram. 
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FIG. 56: Jet-pair invariant mass distribution: the number of events/2 GeV bin after imposing the full set of cuts discussed 
in the text for the benchmark model SPSla'. RH(LH) beam polarization is employed in the top(bottom) panel, assuming an 
integrated luminosity of 250 fb~^ for either polarization. The SM background is shown as the black histogram. 
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FIG. 57: Jet-pair energy distribution: the number of events/2 GeV bin after imposing the full set of cuts discussed in the text 
for the benchmark model SPSla'. RH(LH) beam polarization is employed in the top(bottom) panel, assuming an integrated 
luminosity of 250 fb~^ for either polarization. The SM background is shown as the black histogram. 
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Radiative Chargino Production 

As was pointed out by Gunion and collaborators |3^, 
|37| , Q (see also Riles et a/.||3^), in the case where Am^ 
is in the approximate range 0.1 GeV Am <, 2 GeV, 
the dominant decay mode of the charginos is into soft 
pions plus the LSP (which appears as missing energy), 
as shown in Fig. The dominant SM background to 
this final state is from 77 interactions and has an enor- 
mous event rate. If Am^ < m^r, however, the dominant 
chargino decay is leptonic, xt ~^ ^^^eXij ^iid does not 
pose a problem for detection as this essentially results in a 
charged stable-particle search. In this subsection, we dis- 
cuss the radiative chargino search where the hard photon 
emitted in the process e+e" xtxil is tagged. Unlike 
the other chargino signatures we consider, the strength 
of this signal is not very dependent on the mass splitting 
between the lightest chargino and the LSP neutralino. 

Event Generation 

An immediate issue in performing the search for radia- 
tive chargino production is that we find PYTHIA under- 
estimates the rate and energy distribution of hard pho- 
ton emission from the final state charginos. Thus we use 
CompHEP to generate the e+e~ XiXil well 
as e'^e' — * ^^Xi) events from the explicit (tree-level) 
matrix elements. 

In particular, we find that PYTHIA with defauft ISR 
and FSR options yields a lower cross section for chargino 
production with an associated photon than does Com- 
pHEP. This is illustrated in Fig. ^ where we display the 
cross section for e+e~ X^X 7 in one of the AKTW 
models (labelled as model 13348) as a function of photon 
transverse momentum. For this model (where the 
has a mass of ~ 124 GeV) the cross section computed by 
PYTHIA is about 20% smaller than the CompHEP cross 
section for all values of pT- In examining the MSSM pa- 
rameter space further, we find points where the PYTHIA 
generated cross section can be as low as 50% of that cal- 
culated via CompHEP; for the models considered in this 
paper, the PYTHIA cross section is generally 80 — 90% of 
that from CompHEP. As CompHEP uses an exact (tree- 
level) matrix element calculation, it is presumably more 
accurate. Therefore we use CompHEP when calculat- 
ing the cross sections for radiative chargino production 
in each of our models and to generate the events for this 
process. 

CompHEP does not allow one to set an arbitrary beam 
polarization. However, one may set the electron beam 
polarization to be, for example, purely left-handed, by 
effectively inserting the relevant projection operator into 
the expression for the matrix element. Thus we can cal- 



culate the desired cross sections and generate events for 
each of the two initial helicity states that we consider 
here. For each pure initial helicity state we generate two 
large event files and find the relevant cross sections. We 
then choose the correct number of events for each of our 
two partial electron beam polarizations (80% left-handed 
and 80% right-handed) and pipe this through our analy- 
sis as described in Section . 

CompHEP includes the option of using a beam- 
strahlung spectrum calculated from the beam dimensions 
and the number of particles per bunch. However, for con- 
sistency with the rest of our signal and background, we 
must use the same beam spectrum as described above. 
To implement this spectrum in CompHEP, we read in 
the beamspectrum as generated by GuineaPig ||^ Q . We 
checked that the normalization was correct by compar- 
ing the cross sections for e^e~ — > tt as generated by our 
modified CompHEP code and by PYTHIA. 

In calculating the cross section and generating events 
for the Xi^xr7 final state, we demand that the transverse 
momentum of the photon be greater than 5 GeV. This 
cut is much softer than that we apply in our analysis with 
the detector simulation; we do not wish to eliminate the 
possibility of low pr signal events passing the final pr cut 
due to mismeasurement. We do, however, need to apply 
a cut at this stage for the purpose of regularization. 

Analysis 

In our analysis, we tag on a high-pr photon, produced 
by the signal either off the initial state electron-positron 
pair, or radiated off of one of the charginos. We apply 
the following kinematical cuts as suggested in ||3^, |3^ : 

1. There be exactly one photon in the event with 
Pt > 0.035-^5 and no other charged tracks within 
25 degrees. This isolation cut removes most of the 
two-7 background. 

2. There be no identified {i.e., above 142 mrad) elec- 
trons or muons in the event. Although this cut 
slightly reduces the signal, we find it dramatically 
decreases the background from 77 and 6^7 events. 
For the signal, we see from Fig. ^that the branch- 
ing fraction of charginos to electrons or muons is 
less than 30 — 40% in the relevant Am^ range. 

3. We demand that the number of charged tracks 
be in the range 1 to 11. Note that below 142 
mrad the detector only observes clusters of en- 
ergy, however, we nonetheless treat clusters on the 
same footing as tracks. This cut removes high- 
multiplicity events. In particular, the removal of 
high-multiplicity events restricts this analysis to 
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FIG. 58: Comparison of the cross section for the pair production of the Ughtest chargino with an associated photon as a function 
of the photon transverse momentum as calculated by PYTHIA (red, bottom curve) and by CompHEP (blue, top curve) in one 
of the AKTW models where m = 124 GeV. 



the range of Am^ that we are targeting in this anal- 
ysis. Models with larger values of Am^ generate 
harder partons in the chargino decay that radiate 
more gluons and hence result in more tracks. 

4. We demand that the photon energy and the 
energy of the remaining visible particles satisfy 
-BviB, other particles — < 0.35^/s. This cut furthcr 



mass, which we take to be 160 GeV from the ap- 
proximate 80 GeV lower limit on the chargino mass 
from LEP II ||. 

After applying these cuts we examine the recoil mass 
of the tagged photon. 



M„_, = (1 - 2£;^/^^) . (25) 



reduces the two-7 background. It also serves to re- 
strict this analysis to the relevant range of Am^ , as The dominant remaining SM background again arises 



the amount of visible energy increases with Am^ . 

5. We demand that the ratio of total visible transverse 
momentum to transverse energy satisfy > 0.4 

and that the ratio of total visible transverse mo- 



mentum to total momentum be 



> 0.2. This 



removes most of the hadronic two-7 ^-nd e~^e ini- 
tiated processes. 

6. We require that the recoil mass be Af,ecoii = 
- 2E^l^s) > 160 GeV. This is the recoil 
mass of the tagged photon, which should be at 
least twice the current lower bound on the chargino 



from the reaction e'^e —> I vivid'^ as illustrated in 
Fig. g. 

This analysis is designed to only catch charginos in a 
relatively narrow mass and Ato^ range. If Am-j^ is, e.g., 
greater than ^3 — 4 GeV, then the kinematical prop- 
erties of the average chargino decay may have difficulty 
satisfying our energy, momentum and multiplicity cuts. 
If the mass splitting is too small (less than ~ 0.15 GeV), 
then the chargino will have a long lifetime, as shown in 
Fig. and will decay to at most 2 charged tracks and 
will not pass the above cuts. Furthermore, if the mass 
of the chargino is too close to the beam energy (> 225 
GeV), then not only will the cross section be phase space 
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FIG. 59: Recoil mass of the tagged photon for the SM background to radiative chargino production after the cuts hsted in 
the text have been apphed. This is generated from 250 fb~^ of SM events with 80% right-handed (sohd blue line) and 80% 
left-handed (solid red line) electron beam polarization and unpolarized positron beam at y/s = 500 GeV. The dotted lines 
show the main processes contributing to the background, e"'"e~ ^ l~ viVid''^ , for 80% right-handed (dotted pink line) and 80% 
left-handed electron polarization (dotted green line) . 
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FIG. 60; Chargino lifetime as a function of the chargino-neutralino mass splitting, Am. 
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suppressed, but it will be almost impossible for the sig- 
nal to pass the required photon pr cut. Depending on 
exactly how these kinematic boundaries are chosen, we 
find that only ~ 26 of the 53 AKTW models with kine- 
matically accessible charginos at — 500 GeV have 
these necessary properties. From this analysis we find 
that there are only 14 models which are observable in 
the radiative channel (for either beam polarization) over 
the SM background with a significance 5 > 5. Note that 
although the backgrounds are larger in the LH sample 
than with RH beams, as is usual, the chargino signal in 
this case is far larger (by approximately a factor of 9) 
in the LH sample. This is because for small Am-^ the 
charginos in these models are mostly wino in content. 
This is illustrated in Fig. |6l| which shows our analysis re- 
sults for a number of sample representative models with 
either beam polarization. Model 38239, which is shown 
in this Figure, provides a nice example of a case that is 
missed by this analysis; this model has /^rriji — 0.45 GeV 
while m~+ = 239.75 GeV and thus has little remaining 
phase space to allow for the emission of a hard photon 
with pT > 17.5 GeV. 

The benchmark model SPSla' also leads to a reason- 
able signal excess in this channel as shown in Fig. |6^ . 
However, in this case, this is not the result of a small 
Xi — Xi mass difference, which is 86 GeV in this model, 
but rather the signal in this channel for is a 'fake' in- 
duced by X2 production. As discussed above for some 
of the AKTW models, the observed signal in this case 
is actually a feed down from the production of other, 
perhaps more massive, states in the SUSY spectrum as 
well as from radiative associated X1X2 production. How- 
ever, the xf in SPSla' is already clearly observable in 
the other channels discussed above. We generally find 
that the fake contamination in the radiative channel is 
less than 30% and is quite a bit smaller in many cases. 



77 



M(chi-i-,chi-) 



— L13960.aida 
L16988.aida 

— L24230.aida 

— L27712.aida 
L29820.aida 

— L38239.aida 

— mixLI .aida 

LI 3960.aida 
Entries : 4293 
Mean : 399.84 
Rms : 78.295 
LI 6988.aida 
Entries : 3374 
Mean : 412.67 
87.438 
L24230.aida 
Entries : 4651 
Mean : 397.12 
77.503 
L27712.aida 

Entries : 5545 
Mean : 390.65 
78.233 




260 280 300 320 



M(chi+,chi-) 



— R13960.aida 

— R16988.aida 

— R24230.aida 

— R27712.aida 
R29820.aida 
R38239.aida 

— mixR1.aida 



Ri aaeo.aida 
Entries : 1044 
Mean : 387.53 
Rms : 94.930 




260 280 300 



FIG. 61: Recoil mass distribution for chargino pair production in the AKTW models from the photon tag analysis: number 
of events/2 GeV bin assuming an integrated luminosity of 250 fb"^^ for both LH (top panel) and RH (bottom panel) electron 
beam polarization. 



78 



M(chi+,chi-) 



32 
30 
28 
26 
24 
22 
20 
18 
16 
14 
12 
10 
8 

6- 
4 



MM 




260 280 300 320 340 360 380 400 420 440 



M(chi+,chi-) 



■ L-SPS1aPrimeSB,a 

■ mixLl.aida 



H 1 h 



260 280 300 320 340 360 380 400 420 440 



FIG. 62: Recoil mass distribution for chargino pair production in SPSla' from the photon tag analysis: number of events/2 GeV 
bin assuming an integrated luminosity of 250 fb^^ for both LH (top panel) and RH (bottom panel) electron beam polarization. 
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Very Close Mass Case 

As mentioned above, when Am^ < the decay 
length of the chargino is long compared with the detector 
size, as is shown in Fig. In this case, chargino pro- 
duction may be detected by searching for two massive, 
essentially stable, charged particles that traverse the full 
detector. Seven of the 53 AKTW models with kinemat- 
ically accessible charginos at the 500 GeV ILC fall into 
this category. To perform this search we demand: 

1. There be only 2 massive charged tracks in the event. 

2. There be no tracks, or energy clusters, within 100 
mrad of the beam. 

3. /3 = ^ < 0.93 for both charged tracks. (This value 
is based on constraints from LEPII ||^.) 

2 

4. The energy of the two tracks satisfy '^2 Ei > 

i=l 

0.75Vi. 

The last two cuts remove most of the background from 
the production of muons. After these cuts are imposed, 
the remaining background should be small, aside from 
detector fakes and possible tails from muon production 
due to detector smearing. We then study the /3 = ■§ 
spectrum for both charged tracks and look for excesses 
in the region of low /?. 

In this analysis we search for stable charged tracks in 
the final state whose energy can be reconstructed via 
a dE/dx measurement. /3 = p/E should be signifi- 
cantly less than 1, which would allow us to easily dis- 
tinguish such tracks from those produced by Standard 
Model particles. However, in the current public version 
of org.lcsim dE/dx is not yet implemented pl[ . We 
therefore employ a cheat algorithm, and smear, by a ran- 
dom amount, the energy of all final state tracks which 
we take from the PYTHIA input before detector sim- 
ulation. The width of a random Gaussian fluctuation 
should mimic the resolution obtainable from a more re- 
alistic TOF (Time Of Flight) or dE/dx measurement. 
There is not yet full agreement among the ILC detector 
experts as to the attainable precision which may be pos- 
sible in the determination of /3 ||4C[| , so we perform two 
analyses, one with a 5% and one with a 10% assumed res- 
olution on (3. (Note that an energy smearing of 5(10)% 
translates into a resolution on /3 of roughly 5(10)%.) We 
note that both the ATLAS and CMS detectors have ex- 
cellent /3 resolutions of 5(3)% respectively, and we 
anticipate that any ILC detector should have a compa- 
rable precision as demonstrated in Ref . |Q . 

As shown in Fig. |6^, the background is indeed negli- 
gible for an energy smearing of 5%. However, some SM 
background from Bhabha scattering with missing energy 



due to initial state radiation and beamstrahlung where 
the forward photon is not detected, leaks into the anal- 
ysis when an energy smearing of 10% is assumed. The 
background displayed in Fig. ^ is almost entirely due to 
this process. Nevertheless, we expect the ILC to have an 
energy resolution better than 10%, so this background 
source should not be a problem. 

Out of the 53 AKTW models with kinematically acces- 
sible chargino pairs, only 7 have values of Am^ < and 
have effectively stable charginos as far as collider detector 
measurements are concerned. As in the previous analysis 
for radiative chargino production, the backgrounds are 
similar for both polarizations, however the chargino pro- 
duction cross sections for these models are about a factor 
of 9 larger in the case of LH polarization than in the RH 
case. This is because charginos with small values of Am^ 
are mostly wino in the AKTW models, corresponding to 
large values of the /i parameter. Figure 64 shows these 7 



models for both values of the assumed (3 resolution and 
we see that in either case all of these models are clearly 
visible above background. 

Of course this search will not just find charginos, but 
will also detect any sufficiently long-lived charged parti- 
cles in the kinematically accessible mass range. In fact, 
out of the 28 AKTW models with kinematically acces- 
sible f 's, we find 3 models (labelled as 207, 27285 and 
29334) which have staus with very long lifetimes; in one 
of these cases (27285) the f is the NLSP. Such a situa- 
tion can occur, e.g., in gauge mediation, where the r will 
decay with a very long lifetime by emitting a gravitino 
p3[ . Of course in such a model x? is not the actual LSP; 
thus the cosmology is non-standard, but viable. With 
the above analysis, these models should also lead to ob- 
servable signals in this channel. 

These two possible candidate for stable charged par- 
ticles are easily distinguished by their angular distribu- 
tions, i.e., spin-0 vs. spin-1/2, as well as by their response 
to the various electron beam polarizations as we will see 
below. 

A similar search can be performed in the case of the 
three models with long-lived staus; the results are seen 
in Fig. Recall that in the case of stau pair production 
the cross section is not only controlled by the f mass but 
also by the tlr mixing angle which governs the stau cou- 
pling to the Z boson. Note that the event rates shown 
here are significantly lower than those of the long-lived 
charginos and so the SM background is potentially more 
serious. It is clear, however, that in the case of a 5% 
resolution on j3 the staus in these 3 models will be ob- 
servable; the situation is more difficult to assess by eye 
in the case of only a 10% resolution. A detailed statisti- 
cal study, however, shows that these 3 stau models will 
lead to signals at the level of significance > 5 for both 
choices of the electron beam polarization and for either 
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FIG. 63: (3 = p/E distribution of the SM background after the kinematic cuts listed in the text have been imposed. This is 
generated from 250 fb~^ of SM events generated with 80% right-handed (sohd blue line) and 80% left-handed (solid red line) 
electron beam polarization and unpolarized positron beam at y/s = 500 GeV. In the Figure in the left panel, we assume an 
energy resolution of 5%, the Figure on the right is for an energy resolution of 10%. 



assumed value of the (3 resolution. We observe that stau 
production in these 3 special models with both LH and 
RH beam polarization lead to comparable cross sections. 

We note that there are no long-lived charged particles 
in the case of SPSla'. 



Summary of Chargino Analyses 

Here, we collect and summarize the results of the var- 
ious chargino analyses presented in this section. We re- 
mind the reader that the set of AKTW models contain 
53 models with kinematically accessible charginos at the 
500 GeV ILC. The critical parameter that determines the 
open decay channels for the chargino, and hence governs 
the appropriate search analysis, is Ato^, the mass differ- 
ence between the lightest chargino and the Xi LSP. Of the 
53 models, 7 have values of Am^ small enough to render 
the lightest chargino essentially stable and it traverses 
the full detector before it decays. These models are visi- 
ble in our stable charged particle search. An additional 7 
models have chargino-LSP mass differences in the range 
Am^ < 1 GeV, and result in final states with large values 
of missing energy plus several soft pions. These models 
are targeted by our radiative chargino analysis. Thirty- 
seven of the models have mass differences in the range 



1 < Am-^ < 6 GeV and the charginos decay into off-shell 
W bosons. For this region, we designed a multi-pronged 
search strategy using 11 observables in 3 decay channels 
(missing energy -I- ^fi, 2 jets /i, 4 jets). Lastly, 2 of the 
models have Arrij^ > M\y and in this case the charginos 
decay into on-shell W bosons. We developed only one 
search analysis here, utilizing the 4-jet mode, and found 
this to be a very clean channel for detecting chargino 
production. 

A summary of how many models are visible above 
the SM background for each observable in our chargino 
search analyses is presented in Table IIL Here, we em- 



ploy, as always, our visibility criteria that the signal sig- 
nificance S > 5. We see that the channel with the mixed 
final state of 2-jets + muon + missing energy and the 
missing energy observable in the 4-jet -I- missing channel 
yield the highest number of observable models and thus 
are the best channels for detecting chargino production 
in the randomly generated AKTW models. 

Figure ^ displays the location of each AKTW model 
with Am^ < Mw in the Am^— chargino mass plane. 
The color coding of the model marker indicates whether 
it is observable in any of our analyses for either beam po- 
larization as labelled in the Figure caption. The location 
of the various model points in this plane reveals the kine- 
matic properties targeted by each search technique. For 
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FIG. 64: Velocity (=/3 = p/E) distribution for long-lived charginos assuming an integrated luminosity of 250 fb ^ and LH 
electron beam polarization. The top(bottom) panel corresponds to a resolution of 5(10)% on /3. 
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FIG. 65: Velocity(=/3 = p/E) distribution for long-lived staus assuming an integrated luminosity of 250 fb 
electron beam polarizations as labelled. The top(bottom) panel corresponds to a resolution of 5(10)% on /3. 
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TABLE III: Number of models that are visible above the SM background with a significance <S > 5 in each observable at 
= 500 GeV with 250 fb~^ of integrated luminosity for each electron beam polarization. 



example, the radiative production analysis captures the 
models with low Am^ as it was designed to do, and also 
detects charginos that arc light enough to be produced 
with a hard photon. Here, we see that all 7 of the essen- 
tially stable charginos are captured by our stable charged 
particle search, 3 models arc only observable via radia- 
tive chargino production, 11 models are visible in both 
the radiative production channel and at least one of the 
off-shell W analyses, while 26 models arc detectable in at 
least one of the off-shell W analyses but not in radiative 
production. 4 out of the 53 AKTW models are not ob- 
servable in any of our analysis channels. However, each 
of these 4 models have a lower cross section due to phase 
space suppression. 
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FIG. 66: Distribution of the chargino-LSP mass difference versus the chargino mass for the AKTW models with Am<i < 6 GeV 
for the xt states that are accessible at = 500 GeV. The blue crosses represent models that are observable in our suite of 
analysis channels based on the xf decay via off-shell W bosons. The green crosses correspond to models that are only visible 
in the radiative chargino production analysis channel, while the magenta ones represent models that yield observable signals 
in both the radiative and off-shell W channels. The black crosses are models that are visible in the stable chargino analysis. 
The red points are the 4 models where the xt state is not observable in any of our analysis channels, essentially due to phase 
space restrictions. 
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NEUTRALINO PRODUCTION 

The neutrahno sector of the MSSM is the most com- 
plex as the mass eigenstates are admixtures between 
the Bino, neutral Wino, and two neutral Higgsino weak 
states. Neutralinos, x^, can be pair produced in e'^e~ 
collisions via two distinct mechanisms: s— channel Z 
boson exchange can make a neutral Wino plus a Hig- 
gsino, while i, u— channel selectron exchange can produce 
Binos and neutral Winos in all combinations. These 
mechanisms ensure that all of the ten possible processes 
e+e" XiXj are potentially accessible with rates de- 
pending upon the sparticle masses and the various mix- 
ing angle factors. At a •\/s=500 GeV machine, it is likely 
that only the first one or two of these states will be kine- 
matically accessible and this is indeed the case for the 
wide selection of AKTW models analyzed here as shown 
in Fig. ^ and Table |. If the mass separation between 
Xi and X2 is sufficiently large, then the decay channel 
X2 — * Zxi niay lead to a clear signal if the Z boson is 
not too far off-shell. Unfortunately, in the models exam- 
ined here, if X2 is sufficiently light to be produced it is 
very close in mass to xi and we find that such decays are 
almost impossible to observe. That being the case, we 
only consider Xi pair production with a radiated photon, 
as well as X1X2 associated production in the discussion 
below. Recall that xf is the only accessible MSSM par- 
ticle in many of the AKTW models. 

It is important to consider the weak eigenstate mixture 
of the Xi in our set of models; this is shown in Fig. |6^, 
where it is interesting to observe that the lightest neu- 
trahno is mostly a pure weak eigenstate. 

X2X1 Associated Production 

In order to get a handle on the neutralino sector, it is 
important to consider the associated production of neu- 
tralinos, i.e., e"'"e~ — > X2X11 which proceeds by Z boson 
exchange in the s-channel via the Wino and Higgsino 
content of the Xi 2 and by selectron exchange in the t, u- 
channels via their corresponding Wino and Bino content. 
The cross section for this process is thus sensitive to the 
mixing in the neutralino sector as well as to the masses 
of the exchanged cl.r- Note that if the selectrons are 
heavy, as is the case in many of the AKTW models, then 
the s-channel transition dominates; in this situation the 
associated production process will be suppressed if either 
or both of the X2 or Xi have a large Bino content, which 
as we saw above, is a relatively common occurrence in 
the models considered here. 

At Vs = 500 GeV, 46/242 of the AKTW models have 
the final state X2X1 kinematically accessible. The state 
X2 can decay in several ways depending on the mass spec- 



trum details in the gaugino sector. A mode which is 
always present and yields a relatively clean signature is 
X2 XiZ/H, with the Z/H being on- or off-shell de- 
pending on the X2 ~ Xi mass difference. Certainly, this 
channel will be easier to observe in the on-shell case since 
the invariant mass distribution of the visible particles in 
the final state will be peaked at the Z/H mass. In ei- 
ther case, we consider the decay modes Z^H ^ jj, with 
the jets not flavor tagged, as well as the leptonic modes 
Z —)■ e^e~. In order to access the viability of this 

channel, we examine in Fig. ^the mass splitting between 
the first two neutralino states for the 46 models where 
this production mechanism is kinematically accessible. 
Here, we see that for most models the mass splitting is 
rather small; only 8 of these models have neutralino mass 
splittings larger than Mz and 13 have mass differences 
larger than 20 GeV. It is unlikely that any of the other 
remaining models will produce hard enough jets and/or 
leptons to pass the analysis cuts or be visible above back- 
ground. For the models with the larger mass differences, 
we stress again that their signal rates will be controlled 
by both the selectron masses and the Bino content of 
the two neutralinos. The dominant background we con- 
tend with arise from, e.g., e+e^ ZZ + i^i^, 
76 vW with W jj, as well as 77 t^t" . 

To reduce the SM background for associated neutralino 
production, we demand: 

1. There be precisely one lepton pair (electrons or 
muons) or one jet-pair in the event and no other 
visible particles. 

2. The missing energy satisfy E,-^,^, > 300 GeV. This 
removes the majority of the background arising 
from Z and W boson production. 

3. The transverse momentum for each lepton or jet 
satisfy pr > 0.14y^. This cut removes most of the 
ubiquitous 77 and ej initiated backgrounds. 

4. The angle between the lepton or the jet pair be 
< 95 degrees. This further reduces the background 
from W boson production. 

We then examine the invariant mass spectrum of the 
electron-, muon-, or jet-pair. The remaining background 
after these kinematic cuts are imposed is displayed for 
the /i/i , jj+ missing energy channels in Figs. ^ and |7^; 
we find that the background for the e+e^ final state is 
qualitatively similar to those for muons. 

The signal should have a clear peak in the invari- 
ant mass spectrum that reconstructs to the Z boson al- 
though excesses may also appear elsewhere in the distri- 
bution. Note that jet energy resolution is crucial here 
as some background sources, for example, e~7 — > Vf.du 
and e~7 VeSc, have an invariant mass peak at the 
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FIG. 67: The composition of the lightest neutralino in the 242 AKTW models. 
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FIG. 68: X2 ~ Xi mass difference for the 46 AKTW models which have X2X1 associated production kinematicaUy accessible at 
J~s = 500 GeV. 




FIG. 69: SM background to associated neutralino production for the Mi„„(/i/x) distribution. This is generated for a 250 fb~^ 
sample of SM events with 80% RH(solid blue) or LH(solid red) electron beam polarization at = 500 GeV. The dotted 
green line represents the dominant contribution after the cuts are imposed, 77 i~ ^ which is independent of the beam 

polarization. 
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FIG. 70: SM background to associated neutralino production for the Mi„v{jj) distribution. This is generated for a 250 fb~^ 
sample of SM events with 80% RH(soIid blue) or LH(soIid red) electron beam polarization at -^s — 500 GeV. The dominant 
contribution after the cuts are imposed arise from 76 — > u^qq and are shown by the green(pink) dotted line for LH(RH) beam 
polarization. 



W boson mass. Thus the W and Z boson mass peaks 
must be separable in the 2-jet channel. As is common 
in many of our analyses, the SM background is far lower 
with RH electron beam polarization as this suppresses 
W boson production. Note that we may also have to 
deal with backgrounds arising from other SUSY produc- 
tion processes that can fake the signals from associated 
production. 

Typical results for these analyses are shown in Fig. |^ 
for representative AKTW models. In the case of the di- 
jet analysis, three peaks are observable at the masses of 
the W, Z and 120 GeV Higgs boson. Some models lead 
to small excesses on the W peak while some have ex- 
cesses at the Z; others have excesses at both locations. 
Five models are found to show a signal with a signifi- 
cance > 5 in the dijet channel; all these models have 
small excesses at the Z peak, corresponding to the as- 
sociated production channel under consideration. Simi- 
larly, some models also show some excess at the W peak 
arising from a different X2 decay channel: X2 ~* ^^xf 
with W — > jj, Xi ~^ Xi+ very soft jets. This can hap- 
pen in models with light charginos which have a small 
mass splitting with the LSP. Unfortunately, the rest of 
the AKTW models are unobservable, being buried in the 
dijet case by the enormous W jj peak. There are two 
ways to reduce this background: either decrease the jet 
pair mass resolution to a value below 30% /^/E and/or 



employ positron polarization to reduce the SM rate for 
76 Wv. 

In the dimuon channel, the signal region is seen to have 
very little background, however there are also very few 
signal events. A total of ten AKTW models are found to 
show an excess over background with a significance > 5. 
However, only a few of these excesses can be seen in the Z 
mass region. Unfortunately, all but 2 of these models are 
fakes in the sense that they do not have the X2X1 chan- 
nel kinematically accessible; they do, however, all have 
visible smuons. Some of the signal for models which pop- 
ulate the lower invariant mass region originate from addi- 
tional sources, such as XiXi , or even X3X17 production. 
It would appear from these results that perhaps the cuts 
employed in this analysis are too strong even though the 
signal region is essentially background free. However, we 
find that relaxing the cuts, even just slightly, overwhelms 
the signal region by background. We have not found a set 
of cuts that allows more of the signal to be visible over 
background in this channel. A similar situation happens 
for LH electron beam polarization. While some models 
lead to observable signals in the dimuon channel, they 
are all fakes in the case of either polarization and the 
apparent signal is due to feed down from other SUSY 
sources. In the e^e~ channel, 7 models are observable 
with a significance > 5, however only one of them is not 
a fake as illustrated in Fig. It^. We note that there are 
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FIG. 71: Invariant mass distribution in the dijet (dimuon) channel in the top (bottom) panel from the analysis for associated 
neutralino production for representative AKTW models: events/2 GeV bin assuming RH polarization and 250 fb~^ of integrated 
luminosity. As usual the SM background corresponds to the black histogram. 
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fewer fake signals in the dijet channel. 

Thus at this level of statistics, these jj, iJ.'^n^, and 
e+e" analyses have captured very few of the AKTW 
models where X2X1 is kinematically accessible. 

For the conventional benchmark point SPSla', associ- 
ated neutralino production at y/s = 500 GeV can easily 
proceed as the mass of X2 is only 184 GeV. In this case, 
the X2 — Xi mass splitting is ~ 86 GeV, i.e., < Mz, and 
thus the signal is not observable in the dijet channel due 
to the very large SM W boson background. However, a 
reasonable non-resonant signal excess is observable over 
background in the substantially cleaner dimuon mode. 
This is true for either electron beam polarization, how- 
ever the signal is more strongly observable for the case of 
LH polarization, as can be seen in Fig. |7^. Feed down to 
this final states from the production of heavier chargino 
and neutralino states is also present in this model. 
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FIG. 72: Invariant mass distribution for the dielectron channel in the top (bottom) panel for RH (LH) polarization fr 
the analysis for associated neutralino production for representative AKTW models: events/2 GeV bin assuming 250 fb^^ 
integrated luminosity. As usual the SM background corresponds to the black histogram. 
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FIG. 73: Invariant mass distribution in the dijet (dimuon) channel in the top (bottom) panel from the analysis for associated 
neutralino production for model SPSla': events/2 GeV bin assuming RH polarization and 250 fb~^ of integrated luminosity. 
As usual the SM background corresponds to the black histogram. 
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Radiative Neutralino Production 

In 91 of the 242 AKTW models the neutrahno LSP 
Xi is the only kinematically accessible SUSY particle at 
y/s — 500 GeV. The process e+e" 2xi is, by itself, 
impossible to observe as the final state particles are sta- 
ble, neutral and weakly interacting. The only way to 
render Xi production observable is to tag it by the emis- 
sion of a photon off' the initial state electrons or off the 
intermediate channel selectron; one then looks for an 
excess of events of the form e^e~ — > 7 -I- Emiss- The 
SM background to such a signature can be quite large 
and arises mainly from the reaction e+e^ — > vD^ which 
occurs through intermediate W and Z boson exchanges. 
As will be discussed below, beam polarization can play 
an important role in reducing this dominant component 
of the SM background, as W bosons couple to electrons 
in a purely left-handed manner. 

As noted above, we employ PYTHIA for the generation 
of the signal events However, PYTHIA does not take into 
account photon emission from the virtual t-channel selec- 
tron in neutralino pair production. Without this contri- 
bution to the cross section for e^e~ XiXil^ the signal 
would always be invisible beneath the background pro- 
vided by e^e^ vv^. Thus an accurate modelling of the 
radiative LSP signal at the ILC requires a more sophis- 
ticated approach. We therefore use CompHEP to gener- 
ate the full matrix element for this process, in a manner 
analogous to that described above in Section 5.2.1 for the 
radiative chargino analysis. The CompHEP evaluation of 
the cross section for radiative LSP production uses the 
complete matrix element and can be up to a factor of 
2 larger than that given by PYTHIA and also generally 
yields harder photons. 

We tag on a high-px photon, which is the sole visible fi- 
nal state particle in the process e+e~ 7XiXi- Clearly, 
right-handed electron beam polarization should be effec- 
tive in reducing the background contributions from W 
boson exchange in vv^ production. In fact, after the cuts 
described below are employed, we find that the RH SM 
background event rate is about a factor of 7-8 less than 
that with LH beam polarization. (We note that in the 
case of 100% electron beam polarization, the LH cross 
section is almost 50 times larger than that for the RH 
case.) In contrast, we find that the signal cross sections 
for the AKTW models follow either one of two patterns: 
{i) the LH and RH polarized cross sections are compara- 
ble in magnitude or (m) the RH polarized cross section 
is far larger than that of the LH case. Thus, for either 
of these possibilities, RH electron beam polarization is 
highly favored in order to increase the signal and reduce 
the background. We will thus limit ourselves to this po- 
larization configuration in our analysis below. 

We employ the cuts of Ref . Q , and require: 



1. There be exactly one photon and no other visible 
particle in the event. 

2. The photon transverse energy satisfy = 
EySinOy > 0.03^/s. Here, 9^ is the angle of the 
photon with the electron beam axis. 

3. The photon be present in the angular region 
cos 6*^ < 0.9 

4. The total photon energy satisfy the constraint 
E^ < 0.5y/s — 90 GeV. This removes radiative re- 
turn to the Z-pole. 

We then examine the photon energy distribution and 
look for a signal in excess of the SM background; some 
typical results are presented in Fig. [t^. Unfortunately, 
as can be seen from this Figure, S/B is at best ~ 8 — 9% 
for the models shown here. This remains true for all 
180 AKTW models that have kinematically accessible 
Xi states. In many cases S/B is far below the 1% level. 
However, we find that 17 of the models lead to a signal 
significance S greater than 5. We note that in 4 of these 
17 models, the Xi is the only kinematically accessible 
SUSY particle at ^/s = 500 GeV. Of course, a priori, one 
cannot be certain that the neutralino LSP has been pro- 
duced and discovered, as this final state may receive siz- 
able signal contributions from other SUSY sources such 
as e"*'e^ — > vv*^, with v — > i^Xi- In fact, model 36022, 
shown in the Figure, is an example of one such case. 
This renders it difficult to uniquely identify the signal as 
arising from only the lightest neutralino without further 
analysis. 

Of course, increased luminosity or an adjustment of 
the cuts may make the signal in this channel slightly 
more visible, but what would be more useful, inde- 
pendently of the choice of cuts, would be to include 
positron polarization [ ^Sf . Having such polarization at 
the 30(45, 60)% level would reduce the background by 
roughly ~ 44(60, 73)% in comparison to that with the 
canonical 80% electron beam polarization assumed in our 
analysis. The corresponding increase in the signal in the 
most conservative AKTW model would be 24(36, 48)% 
and thus significant boosts in S/B would result; these 
increases can be somewhat larger depending upon the 
parameter values in a particular model. 

The observation of radiative neutralino LSP pair pro- 
duction is rather straightforward in the case of the famil- 
iar benchmark model SPSla', where the LSP is rather 
light with a mass of only 97.7 GeV. Figure ^ shows 
that the signal is much larger in this case than in any of 
the AKTW models. In fact, an excess in the number of 
events over background can be observed for almost the 
entire range of the photon energy spectrum. However, 
some of this excess may be attributed to radiative sneu- 
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FIG. 74: Photon energy spectra for several representative AKTW models from the photon-tagged Xi P^ir production process. 
Shown is the event rate/2 GeV bin assuming RH electron beam polarization and an integrated luminosity of 250 fb~^. The 
black histogram is the SM background. The top panel shows signal plus background for models with larger event rates whereas 
the bottom panel displays more typical cases with signal and the background now being shown separately. 
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trino pair production which is reasonably significant 
this model [|||. 




FIG. 75: Photon energy spectra for the benchmark model SPSla' from the photon-tagged xi pa^ir production process. Shown 
is the event rate/2 GeV bin assuming RH electron beam polarization and an integrated luminosity of 250 fb"'^. The black 
histogram is the SM background. 
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Particle 


Number Visible 


El 


8/9 


en 


12/15 


flL 


9/9 


p-R 


12/15 


Tl.2 


21/28 




0/11 


Vt 


0/18 


xt 


49/53 




17/180 


X°2 


5/46 



TABLE IV: Number of models, at ^ = 500 GeV, which have 
a given final state particle visible above the SM background 
with a significance 5 > 5 as defined in the text, divided by 
the number of models where the same particle is kinematically 
accessible. 



MODEL VISIBILITY AND COMPARISONS 

As discussed in the Introduction, before addressing the 
issue of model differentiation we first need to determine 
which SUSY particles are visible above the SM back- 
ground in each of the models under study. In particu- 
lar we would like to know how many models contain a 
given SUSY particle that is clearly observable with a sig- 
nificance, S, greater than 5 from our Likelihood analysis 
discussed in Section 3. This information can be obtained 
by combining the individual results found in Sections 4 
through 6; the summary of these analyses in terms of 
SUSY particle discovery is displayed in Table [V. This 



Table shows the number of models with a given SUSY 
particle that we found to yield a visible signal above back- 
ground in our analysis relevant for that particular SUSY 
state, as compared to the number of models where the 
same particle is kinematically accessible. Thus, e.g., the 
Cfj is observable in 12 of the 15 models in which it is 
kinematically accessible. We declare a particle to be vis- 
ible for a given model if it is kinematically accessible and 
a signal with 5 > 5 is observed in the relevant search 
channel. We note that it is possible that some of these 
observable signatures may be due to fakes, i.e., the pro- 
duction of other SUSY states; this is certainly true in, 
e.g., the case of radiative xl production. From this Ta- 
ble we see that for the set of AKTW models the ILC 
does an excellent job at detecting selectrons and smuons 
as well as charginos, however, staus are somewhat more 
problematic, and the neutralino sector appears to be dif- 
ficult. 

We can now combine the results represented in the Ta- 
ble and ask for the total number of models which contain 
visible sparticles with a signal significance greater than 
5. Out of the 85 models which have at least one charged 



SUSY partner kinematically accessible, we find that 78 
have visible sparticle signatures at the ILC. The SUSY 
particles in the other 7 models are not detectable mainly 
due to phase space suppression of the SUSY cross sections 
as discussed in the previous Sections. Of the 96 models 
which have only stable neutral SUSY partners accessible 
{Xi or 4 of them are observable via the photon tag 
recoil analysis. Thus, out of all the models with at least 
one accessible SUSY partner we find that 82/181 lead 
to detectable signals at the ILC. This corresponds to 82 
visible models out of the full set of 242 AKTW models 
(recall 61 of the models have no SUSY partners acces- 
sible at 500 GeV). Surveying all of the models, there is 
a total of 129 charged sparticles which are kinematically 
accessible and we find that 111 are visible in our analyses; 
several more may appear as 'fakes'. 

Using our ILC analyses, we now pairwise compare the 
models that were found to be indistinguishable at the 
LHC by AKTW. We recaU that out of the original 283 
model pairs, 121 were removed from our sample due to 
the PYTHIA feature which shifted the LSP mass, leaving 
us with 162 pairs of models to examine. Interestingly, out 
of these 162 model pairs, 90 have only neutral sparticles 
kinematically accessible in both models. 

In order to compare signals that originate from dif- 
ferent models for the observables described in the previ- 
ous Sections, we perform a analysis of the generated 
histogram distributions for the model pairs. To begin, 
recall that we have generated two complete and statis- 
tically independent background sets, Bl and B2, for all 
of the individual analyses. Taking the set of pure sig- 
nal distributions for the two models we wish to compare. 
Ml and M2, we add each signal distribution to one of 
the corresponding background distributions. This forms 
the combinations R ^ Ml + Bl and S = M2 + B2 for 
each observable. We then perform a x^ analysis of the 
two distributions for each model pair, accounting for the 
fact that the number of events in each sample can be 
different: 



E 



with 



R = 



E 



Ri 



Ri -\- Si 



(26) 



(27) 



Ri and Si denote the number of events in bin i produced 
by the two models (plus background) in each observable 
that we compare. Note that such a x^ test is somewhat 
sensitive to the binning of the data, especially since we 
compare two sets of generated "data" instead of com- 
paring a signal to a theoretical prediction. Note fur- 
ther that the above x^ prescription relatively normalizes 
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the two distributions so that wc only compare shapes at 
this point. We then add an additional term to the 
which accounts for the total number of events in both 
histograms and allows for an 1% systematic error in the 
relative normalizations due to luminosity and cross sec- 
tion normalization uncertainties. 

Wc then compute the distributions for each of the 
model pairs, for each the following observables, taken one 
at a time, which were obtained after applying the analysis 
cuts described in the Sections above: 

• Selectron analysis: Ef,+ „ g- and VtV\^- 

• Smuon analysis: ^- and P^yis ■ 

• Stau analysis: Et and VtVi^- employ the r 
identification procedure described in Section , with 
and without the inclusion of electrons in the final 
state in order to remove the background from beam 
remnants. 

• Sneutrino analysis: missing energy for two chan- 
nels, 4-jet plus lepton-pair and 6-jets. Each channel 
is analyzed with two different assumed minimum 
values of the LSP mass. 

• Chargino non-close mass case: For the case of on- 
shell W boson production in chargino decays, we 
examine the E j^tpair spectrum. For the case where 
the charginos decay into off-shell W bosons, we ex- 
amined three decay channels. Our observable for 
the fully Icptonic channel is ^- . In the fully 
hadronic channel we analyze the missing energy dis- 
tribution for the 4-jet final state and the invariant 
mass spectrum of the two jet-pairs. For these two 
distributions we perform analyses with two differ- 
ent assumed values of the Xi mass and also with 
an additional cut on Vty\s,- semi-leptonic 
channel with the jet-pair +//,+ missing energy final 
state, wc examine Ej^ as well as the invariant mass 
of the jet-pair. In this case, we again employ anal- 
yses with two different assumed values of the LSP 
mass. 

• Chargino radiative production: the recoil mass 
Afpecoii of the tagged hard photon. 

• Chargino very close mass case: /? = of the two 
massive tracks in the event, assuming an energy 
smearing//? resolution of 5% and 10%. 

• X2X1 associated production analysis: the invariant 
mass of electron, muon and jet-pairs. 

• Radiative Xi analysis: E^. 



Before examining the results of our x^ model compar- 
ison, we first check our procedure by comparing the two 
pure background samples to verify that, though inde- 
pendent, they arc not statistically distinguishable. Wc 
do indeed find this to be the case for every observable in 
each analysis. Next, we examine each observable listed 
above and determine whether the comparison probability 
(as given by the value of the x^ and number of degrees of 
freedom) shows a difference at the 5(or 3)(j level for that 
specific distribution. Note that wc perform this compar- 
ison separately for each electron beam polarization since 
we have distributions for each polarization configuration. 
If there is a 5(or 3)cr level difference in at least one dis- 
tribution then we claim that the two models are distin- 
guishable at that level of confidence. In fact, a number of 
observables are distinguishable at this level in many mod- 
els. We can improve this procedure by summing over the 
various observables in the x^ computation, taking only 
one distribution from the different analyses to ensure that 
we do not introduce any effects from correlations. Note 
that in the case where we employ only a single observ- 
able in this comparison there are no issues of statistical 
independence in contrast to when several distributions 
are combined. 

When performing these comparisons, wc find that 
many model pairs are not distinguishable. This happens, 
e.g., in all of the 90 cases where we compare models which 
have only neutral particles kinematically accessible. The 
few models where we observe excess photons in radiative 
x5 production are already differentiated by other analy- 
ses involving charged sparticlc production, and thus these 
cases do not help with distinguishing pairs of models con- 
taining only kinematically accessible neutral sparticles. 
This implies that wc should concentrate on the 72 po- 
tentially distinguishable model pairs where at least one 
member of the pair has at least one kinematically acces- 
sible charged sparticlc. In this case, we find that a large 
number of model pairs are distinguishable at the 5a level 
in several different analyses. 

Based on the criteria above, using our results from the 
single observable comparison procedure described above, 
we find that 55(63)/72 model pairs where at least one 
model has kinematically accessible charged sparticles are 
distinguishable at the 5(3)cr level. These results are based 
solely on single final state comparisons between models. 
Making use of the combined observable comparison pro- 
cedure described above, we find instead that 57(63)/72 
pair of models are distinguishable, which is only a slight 
improvement. 

The model pairs that are found to be indistinguishable 
fall into two broad classes: (i) those where one model 
in the pair has only a kinematically accessible neutral 
sparticle, e.g., the LSP, which is not visible above back- 
ground and the other model contains a sole accessible 
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charged sparticlc which is also difficult to observe. There 
are 7 model pairs in this category that cannot be dif- 
ferentiated at the 5cr level. Examples are models with 
a heavy selection and smuon which have kinematically 
suppressed cross sections and models which only contain 
f states that are also difficult to observe due to small 
production cross sections, (m) The second class consists 
of 8 model pairs where each model in the pair contains 
a single kinematically accessible charged sparticle. In 7 
of the 8 model pairs, it is the lightest chargino state that 
is produced and found to lead to an indistinguishable 
signature. 



CONCLUSIONS 

In this paper, we have performed a systematic and de- 
tailed analysis of the capabilities of the 500 GeV ILC 
(with an integrated luminosity of 500 fb^^ with 80% elec- 
tron beam polarization) to explore the nature of a large 
number, 242, of scenarios within the MSSM. The goal 
of this project was to study general model-independent 
MSSM signatures and then to determine whether the ILC 
could differentiate between 162 pairs of these models, i.e., 
MSSM parameter space points, which were found to be 
impossible to distinguish at the LHC. To do this we first 
had to address the issues of kinematic accessibility of the 
SUSY states, as well as the experimental observability of 
the corresponding sparticle production over the SM back- 
ground. Given the present controversy about the ILC in 
the world-wide community, we feel that such a detailed 
study of the MSSM at the ILC is critical at the present 
time. 

In order to accomplish this task, we employed a com- 
plete set of full matrix element SM backgrounds for 
all 2 2, 4 and 6 parton final states initiated by 
e~^e',je^ and 77 as calculated by T. Barklow using 
WHIZARD/O'MEGA. We made use of both PYTHIA 
and CompHEP for generating the SUSY model signal 
events and employed a realistic beamspectrum generated 
with WHIZARD/GuineaPig. Additionally, we included 
the effects of the SiD detector by implementing a ver- 
sion of the org.lcsim fast simulation. In this analysis 
we assumed that the integrated luminosity was equally 
split between two distinct samples with 80% LH or RH 
electron beam polarization. Analyses were performed on 
many different SUSY channels simultaneously in order to 
probe the charged slepton, sneutrino, lightest chargino, 
LSP and X2 ~ Xi sectors. A universal set of cuts for all 
models was developed. This is the first realistic, model- 
independent analysis of its kind of SUSY signatures at 
the ILC 

Out of the original 242 models only 85 led to the exis- 
tence of a kinematically accessible charged SUSY partner 



at 500 GeV. The remaining models cither had no SUSY 
particles kinematically accessible (61) or only the light- 
est neutralino and/or sneutrino accessible (96). Using log 
likelihood techniques, we found that 78/85 models with 
a charged SUSY partner as well as 4 additional models 
which only had neutral particle states accessible were vis- 
ible above SM background in our analyses. Thus, a total 
of 82/161 models with accessible particles were found to 
be observable at the 500 GeV ILC. In performing our 
analysis, beam polarization was essential in reducing the 
SM background and allowed for distinguishing sparticle 
states in many cases. Some models contained charged 
states that were found not to be visible generally as a 
result of suppressed cross sections due to phase space 
availability. Of the 72 pairs of models where at least 
one member of the pair contains one or more accessible 
charged SUSY partner, our analysis found that 57(63) 
of the pairs could be distinguished at the level of 5(3)c7. 
These results are generally quite surprising. Perhaps this 
is because the "standard lore" about SUSY at the ILC 
derives largely from studies of a handful of very special 
points in MSSM parameter space, whose signatures are 
particularly easy to see. We, on the other hand, have 
studied a wide variety of MSSM parameter points. 

From this analysis it is clear that the ILC with the SiD 
detector does a reasonably good job at observing charged 
sparticles which are kinematically accessible and distin- 
guishing models containing such particles. The major 
weakness, beyond the restricted kinematic reach, is in 
the neutral sparticle sector. This problem may be re- 
solved by employing positron beam polarization as well 
as even more sophisticated analyses. 

For the future we plan to extend this analysis to the 
case of 1 TcV centcr-of-mass energy, and include a stiidy 
of the influence of positron polarization as well as a num- 
ber of detector issues. 
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